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Introduction 



Almost sixty years past the discovery of the neutrino [1], the study of its physics is 
still one of the outstanding and most exciting challenges to the scientific community. 
Thanks to the striking progress of particle physics in the last decades, the neutrino 
interaction mechanisms in vacuum are now well-established. On the other hand, 
little is known about its nature (whether it is a Dirac or Majorana particle), its mass 
or the existence of sterile neutrino. Many constraints on its microscopic features, 
most notably the mass, can be inferred from astrophysical data. Recently, neutrinos 
have been also used as astronomical probes, as they can carry unaltered information 
over cosmological distances [2]. 

Although several questions about neutrinos have been answered in the past few 
years, a number of issues are still obscure. In this context, it is worth mentioning 
the over one hundred papers that appeared within a couple of days after the abrupt 
announcement of the observation of superluminal neutrinos. This claim was later 
withdrawn, but in those days many new intriguing models and theoretical tools were 
developed. The neutrino "gold-rush" seems to be never ending. 

The behavior of neutrinos in matter is far from being well understood, and should 
be regarded as a largely independent challenge, involving issues connected with both 
neutrino physics itself and the physics of the system neutrinos are interacting with. 

Neutrino interactions with matter strongly depend on energy. In this Thesis 
we restrict our focus to low-energy neutrinos, whose interactions with many-body 
systems, such as nuclei or, more generally, nuclear matter, turn out to be critical 
to the description of a variety of different scenarios, from supernovae explosion to 
neutron-star cooling. 

The supernova is one of the last stages of star evolution [3] : when the ignition 
of the natural elements is no longer exothermic, the star starts contracting under 
the gravitational pressure. If the star is very heavy, its mass exceeding 4M0, where 
Mq denotes the solar mass, the electron degeneracy pressure is not able to balance 
gravity, and the core collapses, reaching the typical density of a nucleus. Now the 
core reacts elastically to a further compression, leading to a bounced shock wave that 
ejects the mantle of the star. When the expanding ejecta become gravitationally 
decoupled from the stellar residue, the unshocked core evolves into a proto-neutron 
star. For several seconds it is still optically thick, and neutrinos act as a very efficient 
mediator to transfer energy from the core to the colder external region^. As a 
consequence, the temperature rises and the star gets smaller, reaching a radius of 



^This high efficiency is reminescent of the fact that the impact of an elementary particle on the 
stellar structure is maximized when its mean free path is of the order of the size of the system. 
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about 10 — 15 km. After ^ 50 s, the neutrino mean-free-path exceeds the radius of 
the star, and neutrino flow turns into in a coohng mechanism [4]. It can remain the 
dominant coohng mechanism for a period of time that varies from weeks to several 
thousand years, depending on the stiffness of the equation of state, driving in turn 
the onset of the direct Urea processes. 

Experimentally, the importance of neutrino-driven mechanisms has been shown 
by the analysis of the signal from the core collapse supernova SN1987A, carried out 
by the KAMIOKANDE [5] and IBM [6] collaborations. Most of the emitted energy 
turns out to be carried by neutrinos, the relevant figures being [7] : 

Egrav ^ 10^^ erg , ^ 2.7 x 10^^ erg , 

Erad ^ 8 X 10^^ erg Ek^n ^ 10^^ erg . 

From a quantitative point of view, the understanding of neutrino emission rates 
and propagation in dense matter is regarded as one of the most critical issues, 
particularly in view of the large effort being made to carry out realistic large scale 
simulations of supernovae and neutron star evolution [83, 84]. The systematic 
uncertainty associated with these simulations depends heavily on the values of the 
neutrino- nucleon and neutrino- nucleus cross sections used as input. Many results 
discussed in the literature have been obtained using somewhat oversimplified models 
of nuclear dynamics, although this source of uncertainty has been recently reduced 
adopting more realistic models. 

In this Thesis we explore the effects of nuclear correlations on neutrino interactions 
with the nuclear medium. We will discuss the main coupling mechanisms at low- 
momentum transfer and the role of many-body excitations. In this regime, collective 
modes turn out to be important, and must be taken into account. As a final outcome, 
a quantitative estimate of the neutrino mean free path in dense neutron matter will 
be provided. 

Nuclear interactions are described within the ah initio approach based on realistic 
hamiltonians and the Correlated Basis Function (CBF) formalism. Within this 
scheme, nucleon- nucleon correlations are included in the nuclear wave functions, and 
one can define an effective potential, modified by medium effects, suitable for use in 
standard perturbation theory. Unlike the effective potentials derived within mean 
field approaches, the CBF effective potential reduces to the hare potential in the 
limit of vanishing density, and is therefore strongly constrained by the experimental 
information on the two-nucleon system. 

Many-body effects at low momentum transfer are treated within Landau theory 
of normal Fermi liquids, suitable to describe both coherent and incoherent excitations 
on the same footing. The Landau parameters are derived from matrix elements of 
the CBF effective potential, and turn out to provide a description of static properties 
of neutron matter in very good agreement with that obtained from CBF calculations. 

The Thesis is structured as follows: 

• Chapter 1 outlines the main features of nuclear dynamics and the many-body 
approaches commonly employed to study nuclear matter; 
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• Chapter 2 is devoted to the discussion of the approach based on the CBF 
formahsm and the cluster expansion technique, as well as the derivation of the 
effective potential; 

• In Chapter 3 we briefly review Landau theory of normal Fermi liquids; 

• Chapter 4 reports the Landau parameters obtained from the effective potential 
derived in Chapter 2. The static properties of matter obtained from Landau 
theory and the CBF approach are also compared; 

• In Chapter 5 we discuss the dynamic form factors of neutron matter, entering 
the definition of the neutrino cross section in matter; 

• Chapter 6 is devoted to the discussion of the mean free path of non degenerate 
neutrinos, including its dependence on matter density and temperature; 

• As the treatment of three-nucleon interactions is one of the critical issues in 
the definition of the effective interaction, in Chapter 7, we report the results of 
a recent study of their effects on the binding energies of the oxygen isotopes, 
carried out within a novel approach based on the Green function formalism. 
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Chapter 1 

Overview of nuclear many-body 
theory 



The main chahenge of any many-body approach is to provide an ab initio description 
of the system under study. Starting from the "elementary" degrees of freedom 
and including few-body interactions, one should be able, at least in principle, to 
reach a full understanding of the physical behavior. So forth, the system should 
be modeled without the need of any approximations or the need to include new 
phenomenological parameters. However, this lucky chance is somewhat exceptional 
in physics. Even in vacuum, as in particle physics, this situation occurs at some 
energy scales only, depending on the interaction we are dealing with. Focusing on 
nuclear interactions, the celebrated Quantum Chromodynamics (QCD), universally 
accepted as the fundamental theory of strong interaction, provides an excellent 
scheme to deal with high-momentum transfer processes. On the other hand, it is 
clearly not the right approach to describe the low-energy region, where systems are 
highly non-perturbative and any perturbative expansion is designed to fail. 

In this case, a coarse- graining method is absolutely necessary. One usually 
averages on fast variables (like quarks), and defines new variables that are linked to 
a much larger time scale. In nuclear physics this coarse-grained (averaged) variables 
are already provided by nature: the nucleons. They still represent the actual basis 
for any realistic low-energy calculation, classification and prediction of new bound 
nuclei, see Segre chart of fig. (1.1). Several, quite important features of nuclear 
dynamics can be understood and constrained thanks to description of nuclei in 
terms of nucleons: binding energy, symmetries, saturation properties are only few 
of these. Moreover, thanks to the large amount of data available, the equilibrium 
properties of finite nuclei are well established, and consequently the properties of 
infinite system can be extrapolated with little error. To see how this can be possible, 
let us consider the celebrated semi- empirical mass formula, yielding the binding 
energy as a function of the total number of nucleons A and the number of protons 

z, 

B{Z,A) Ir , ,2/3 {A-2Zf ^ 1 1 

___ ^ _ -a,Al - ae-^ - + Aap-^J . (1.1) 

The first term in square brackets, proportional to A, is called the volume term 
and corresponds to the bulk energy. The second term, proportional to the squared 
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Figure 1.1. Table of Nuclides, also referred to as "Segre Chart". The figure shows only 
part of the table, hilighting the different properties of the included nuclei [10] 



nuclear radius, is associated with the surface energy, while the third accounts for the 
Coulomb repulsion between Z protons uniformly distributed within the sphere of 
radius R. The fourth term, that goes under the name of symmetry energy is required 
to describe the experimental observation that stable nuclei tend to have the same 
number of neutrons and protons. Moreover, even-even nuclei (i.e. nuclei having 
even Z and even A — Z) tend to be more stable than even-odd or odd-odd ones. 
This property is taken into account by the last term in the above equation where 
A = —1, 0, 1 for even-even, even-odd, odd-odd nuclei respectively. In fig. (1.2) the 
value oi B/A for stable nuclei is shown together with the empirical values of the 
coefficients. By inspecting fig. 1.2, we can readily infer some important pieces of 
information. The nuclear interaction length is roughly the same as the radius of light 
nuclei, and smaller than the radius of heavy nuclei, otherwise the binding energy 
would increase with increasing A. The properties of infinite matter can be explored 
in the limit A ^ oc, yielding 

E/A^aY ^ 15.7 MeV for A^2Z (symmetric matter, or SNM) , 
E/A^aY ^ -19.37 MeV for Z = (pure neutron matter, or PNM) . 

The above values indicate that, while SNM is stable, PNM is not, since the system 
unbound [E/A < 0). Note that the fact that pure neutron matter alone is unstable 
does not prevent the occurrence of PNM in a broad region in the interior of neutron 
stars. In an astrophysical object of mass comparable to the solar mass and radius 
^10 km has the gravitational pressure is strong enough confine PNM within the 
core. 

Many other features of the nuclear interactions can be figured out studying the 
static properties of finite nuclei and extrapolating to infinite matter. The most 
important are: 

Saturation of charge densities In fig. (1.3) the charge density distributions of 
different nuclei are plotted as a function of the radial coordinate. It clearly 
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Figure 1.2. Upper panel: Binding energy per nucleon for stable nuclei according to eq. 
(1.1), with ay = 15.67 MeV,as = 17.23 MeV, = 0.714 MeV, aA = 93.15 MeV and 
ttp = 11.2 MeV. Lower panel: the solid line shows the volume contribution to the 
binding energy per nucleon, whereas the surface, Coulumb and symmetry contribution 
are represented by diamonds, square and crosses, respectively. 




Figure 1.3. Charge distribution as function of the radial distance, obtained from elastic 
electron scattering on nuclei. 
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1. Overview of nuclear many-body theory 



appears that the density in the nucleus interior is nearly constant and inde- 
pendent of the mass number A. This means that nucleons cannot be packed 
too tightly, implying in turn that at short distance the nucleon-nucleon (NN) 
force must be repulsive. Assuming that the interaction can be described by a 
non relativistic potential v depending on the interparticle distance r we can 
write: 

v{r) > , |r| < Tc , 
where Tc is the radius of repulsive core. 

Behavior of the binding energy per nucleon The nuclear binding energy per 
nucleon is roughly the same for all nuclei with A > 12. This observation 
suggests that NN interaction has a finite range ro, i.e. that 



v{t) , |r| > ro . 



Properties of mirror nuclei The spectra of the so-called mirror nuclei, i.e. pairs 
of nuclei having the same A and charges differing by one unit (implying that 
the number of protons in a nucleus is the same as the number of neutrons in 
its mirror companion), e.g. ^^N (A=15, Z=7) and g^O (A=15, Z=8), exhibit 
striking similarities. The energy levels with the same parity and angular 
momentum are the same up to small electromagnetic corrections, showing that 
protons and neutrons have similar nuclear interactions, i.e. that nuclear forces 
are charge symmetric. 

The last property is a manifestation of the "isotopic invariance", associated with the 
approximated isospin symmetry of strong interactions^. This mean that the proton 
and the neutron can be viewed as different states of the same particle, described by 
the isodoublet field 




The isospin formalism allows one to describe free protons and neutrons using the 
lagrangian density 

invariant under the SU (2) global phase transformation 



where the aj are constants and the rj denote Pauli matrices in acting in isospin 
space. The proton and the neutron correspond to isospin projections +1/2 and 
— 1/2, respectively. Proton-proton and neutron-neutron pairs always have total 
isospin T = 1, whereas a proton-neutron pair may have either T = or T = 1. The 

^The isospin symmetry would be exact if the proton and neutron masses were equal. It is 
explicitly broken by small electromagnetic corrections and by the mass difference between the u 
and d quark. 



1.1 The two-nucleon system 
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two-nucleon isospin states |T, Mt) can be summarized as follows: 

|i,i> 

|i,o) 

|0,0) 

Isospin invariance implies that the interaction between two nucleons separated 
by a distance r = |ri — r2| and having total spin S depends on their total isospin T 
but not on its projection Mt- For example, the potential v{v) acting between two 
protons with spins coupled to S=0, is the same as the potential acting between a 
proton and a neutron with spins and isospins coupled to 5 = and T — 1. 

1.1 The two-nucleon system 

The two-body interaction is the basic ingredient to perform theoretical calculations 
in both finite nuclei and infinite nuclear matter. However, the derivation of this 
interaction from the fundamental theory, QCD, is still out of reach of the available 
theoretical approaches. As pointed out in the previous Section, the behavior of the 
NN force in the r ^ 0, oc limit can be infrrred and a qualitative trend can be drawn. 
This is of course not enough, since the same static properties can be predicted from 
different models of the microscopic dynamics. 

Several different approaches, described in the literature, have been employed to 
overcome this problem. The left panel of fig. (1.4) shows the preliminary results 
of the ongoing effort aimed at obtaining the NN potential within the framework 
of lattice QCD recently reported in Refs. [11, 12]. The radial dependence of the 
potentials in the ^S^o and ^S^i channels have been obtained setting the pion mass 
to ttIt^ ^ 529 MeV. It is apparent that the expected shape, exhibiting a repulsive 
core followed by an attractive region, is reproduced at qualitative level. On the 
other hand, comparison with the corresponding phenomenological potentials (to be 
discussed below) displayed in the the right panel, clearly shows that the attractive 
component is far too weak. Note that the NN potential reduces to the Yukawa 
one-pion exchange interaction at long distances, while the short-range repulsive 
core is to be ascribed to heavy-meson exchange or to more complicated mechanisms 
involving nucleon constituents and the intermediate range attaction is largely arising 
from two-pion exchange processes. Although the progress of QCD calculations has 
been impressive, the presently available potentials, as honestly stated by the authors 
of Refs. [11, 12], cannot be regarded as realistic. Use of these interactions would in 
fact lead to predict that a number of nuclei, including the deuterium, be unbound. 

Chiral perturbation theory [13] has been widely employed to obtain a 

theoretically sound model of the NN interaction. This approach, originally proposed 
by Weinberg, exploits the Goldstone boson nature of the pion, arising from the 
breaking of chiral symmetry [14]^. As a consequence, the interactions of low energy 



|n, n 
1 



V2 



(|p,n) - |n,p)) . 



^The pion must be actually regarded as a quasi Goldstone boson, as chiral symmetry, besides 
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1. Overview of nuclear many-body theory 



pions are weak, and can be treated in perturbation theory. xPT has been first 
successfully applied to the pion-pion and pion-nucleon systems. Many pion systems 
can be still be described within the perturbative approach, using the formalism of Ref. 
[15]. The extension to the nucleon-nucleon sector is conceptually more problematic, 
due to the nature of the interaction that may lead to bound states. However, 
Weinberg has shown that the non perturbative contributions to the NN scattering 
amplitude are associated with diagrams involving purely nucleonic intermediate 
states, that can be summed up to all orders in perturbation theory solving a 
Lippmann-Schwinger equation. 

Up to now, well-tested two-body potential based on x^T have been obtained 
at N^LO (Next-to-Next-to-Next Leading Order). They involve a relatively small 
number of free parameters, to be determined phenomenologically . 




r[fm] 0.5 1 1.5 2 2.5 



Figure 1.4. Left panel: Central componet of the NN potential obtained from lattice 
QCD [11, 12] setting the pion mass to 529 MeV. It is apparent that, while the expected 
shape, exhibiting a repulsive core followed by an attractive region, is achieved, the 
attractive component is much weaker than that predicted by the phenomenological 
potential showed in the {righ pannel) 



Phenomenological approach 

A different approach to attack the problem of the determination of the NN potential 
relies on phenomenology only. For example, the observation that the only observed 
NN bound state, the nucleus of deuterium, consists of one proton and one neutron 
coupled to total spin and isospin 5 = l,r = 0, respsctively, is a clear indication of 
the spin-isospin dependence of nuclear forces. Furthermore the deuteron exhibits a 
non vanishing electric quadrupole moment, implying that its charge distribution is 
not spherically symmetric. Hence, the NN potential in non-central. 

Besides the properties of the two-nucleon bound state, the large data base of 
phase shifts measured in NN scattering experiments (the Nijmegen data base[16] 
includes ^ 4000 data points, corresponding to the energies up to the pion threshold 



being spontaneously broken, is also explicitly broken by the mass difference between the u and d 
quark. 



1.1 The two-nucleon system 
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{350MeV) provides valuable additional information on the nature of NN forces. 
Phenomenological potentials are usually written as a sum of two terms: 

where v^^ takes into account the long range interaction modeled by the one-pion 
exchange (OPE) potential, while vr describes the interaction at both intermediate 
and short range. Let us consider these two terms in more details, 




400 



Figure 1.5. Left panel: Feynman diagram depicting the one-pion-exchange process, driving 
the NN interaction at large distance. panel Comparison between ^Sq and 
phase shifts resulting from the Nijmegen analysis [16] (open circles) and the predictions 
of the Argonne vi^ and OPE potentials. 



Vtj: In the literature, the one-pion exchange potential has been derived using different 
choices of the tvNN coupling. For example one can write the interaction 
lagrangian in the form 

where qa is the axial coupling constant, that can be extracted from neutron 
beta-decay, and F^^ is the pion decay constant. In the above equation, the 
75 accounts for the fact that the pion field, represented by the T = 1 isospin 
triplet TV = (tt^, tt", tt^) with spin-parity 0~, transform like pseudoscalar under 
Lorentz transformations. This interaction term, involving a pseudovector 
nucleon current, is usually employed in the modern versions of xP^- However, 
it totally is equivalent to the widely used pseudoscalar tt NN interaction term, 
that can be obtained through an integration by parts 

. 2mNgA ^ .J 




Note that this equivalence hold true because the interacting nucleons are on 
the mass shell. The tvNN coupling constant, resulting from the Goldberger- 
Treiman[17] relation is 



g^NN 2mNgA/F^ ^ 14 
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The amphtude of the OPE process, described by the Feynman diagram of fig. 
(1.5) reads 

^s' (^2)75^52 (^2)^.^ bi)75^.i (pi) . . 
Mf^ = -ig^NN k^-ml ' 

where k — p'^—pi — p2 —p'2 •> — k^k^ — k'^ — |kp , Us{p) is the Dirac spinor 
associated with an on-shell nucleon, and 

where 77^ is the Pauh spinor describing the isospin state of the i-th nucleon. In 
non-relativistic hmit, the above amphtude provides the expression of the NN 
interaction, the expression of which in coordinate space is 

2 pHlTrr 

^^ = '-i^ir^■T2){cT,.V){a2■V) — 



- in • T2 ) 



(0-i-0-2) + 5i2(l + ^ + ^ 



e^ 47r , , 



(47r)2 4:171% 3 \ ' \ X X" 

where x — rriT^r. Note that the operator 

5^12 3(0-1 • f)(o-2 • f) - (o-i • 0-2), 

responsible of the non centrality of the OPE potential, is reminiscent of the 
tensor operator describing the non-central interaction between two magnetic 
dipoles. 

The OPE potential provides an accurate description of the long range part 
(|r| > 1.5 fm) of the NN interaction, as shown by the very good fit of the 
NN scattering phase shifts in states of high angular momentum. In these 
states, due to the strong centrifugal barrier, the probability of finding the 
two nucleons at small relative distances becomes in fact negligibly small. The 
notation v^^ refers to OPE potential v^^ stripped of its zero-range component. 

vr At medium and short range many other complicated mechanisms, involving 
the exchange of two or more pions or heavier mesons (like the p and cj, with 
masses rup — 770 MeV and = 782 MeV, respectively) contribute to the 
NN interaction processes. Moreover, when the relative distance becomes very 
small (|r| < 0.5 fm) the composite nature of nucleons is expected to play a 
crucial role. 



The potential obtained from the sum of v^^ and vr is usually written in the form 

ST 

where the indices S and T denote the total spin and isospin of the interacting 
pair, respectively, while P and 11 are the spin and isopsin projection operators, 
whose definition and properties are given in Appendix A. The functions VTs{T^ij) ^^d 



1.2 Three- nucleon interactions 
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vtriTij) describe the radial dependence of the interaction in the different spin-isopsin 
channels, and reduce to the corresponding components of the one-pion-exchange 
potential at large r^j. Their shapes are adjusted in such a way as to reproduce the 
available NN data (deuteron binding energy, charge radius and quadrupole moment 
and the NN scattering data). An altrernative representation of NN potential, based 
on the set of six operators, 

or/' = [X, (CTI • CT2), 5l2] ® [X, (n . T2)] , 

is given by, 

% = E^"(n.M^- (1-2) 

The above potential provide a fairly good description of deuteron properties and the 
S — wave scattering phase shift. However, to achieve a good description of the P 
wave one has to include two additional components, associated with the momentum 
dependent operators defined as 

Olf^^ = L • S [X, (n • T2)] , 

where L is the orbital angular momentum. The full expression of the potential, 
yielding the best available fits of NN scattering data, with x^/datum^ 1, includes 
ten additional operator more, bringing their total number to eighteen 

Qn=9,...M ^ jl2, . a2). (L • Sf] ® [I, (n . T2)] , 

(^n=15,...,18 ^ . ^^^^ ^ ^j.^^^ (^^^ ^ ^ 

where 

Ti2 = 3(ti • r)(T2 • r) - (ti • T2) . 

The 0-=-'-'^- take care of small charge symmetry breaking effects, due to the 
different masses and coupling constants of the charged and neutral pions. 

The calculations discussed in this Thesis are based on a widely employed potential 
model, obtained within the phenomenological approach outlined in this Section, 
referred to as Argonne vis [18]. As an example of the quality of the phase shifts 
obtained from the Argonne vis potential, in fig. (1.5) we show the results for the ^Sq 
and partial waves, compared with the predictions of the OPE model (OPEP). 
We will adopt a static version of the full Argonne vis^ including only the first six 
operators. Use of this simplified potential appears to be justified in the context of 
our study of the dynamic form factors of neutron matter in low-momentum transfer 
regime. 

1.2 Three- nucleon interactions 

It has been long realized that three-nucleon interactions play a critical role in 
determining the properties of both finite nuclei and infinite nuclear matter. The 
exact solution of the Scrodinger equation for the ground state of the three-nnucleon 
system with the Argonne vis potential yields a binding energy per nucleon £"0 = 7.6 
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MeV, to be compared to the experimental value E'exp = 8-48 MeV, while accurate 
calculations of the density dependence of the energy per particle of symmetric nuclear 
matter carried out with any phenomenological NN potential fail to reproduce the 
empirical saturation properties. In order to bring theoretical results into agreement 
with the data one has to add to the nuclear hamiltonian a contribution involving 
the potential Vijk^ describing three- nucleon interactions. 

The theoretical description of three-nucleon potential Vijk has been first discussed 
in the pioneering work of Fujita and Miyazawa [21]. They argued that the main 
contribution comes from a two-pion exchange process in which a nucleon-nucleon 
interaction leads to the excitation of one of the participating particles to a A 
{M/^ = 1232 MeV) resonance, which then decays due to the interaction with a third 
nucleon (see fig. (1.6). 

The recent models of the three-nucleon potentials are usually written in the form 
written in the form 

Vijk = y^k + ^ij^k 5 

where the first contribution is the attractive Fujita-Miyazawa term of fig. (1.6, 
while V^j^ a purely phenomenological repulsive term. The parameters entering the 
definition of three-body potential are adjusted in such a way to reproduce ground 
state energies of and ^He and the saturation properties of symmetric nuclear 
matter [22]. 

Three-nucleon interactions deeply affect both the binding energy and several 
spectroscopic features of nuclei, mainly close to the neutron drip line. In this case, 
the last occupied nuclear orbital are in fact very close to the continuum spectrum, 
and the probability of many-body interactions is larger. As an example, in the 
left panel of figure (1.7) we show the binding energy of several isotopes of Oxygen, 
obtained using the Self Consistent Green Function (SCGF) formalism [19]. The 









A 







Figure 1.6. Diagrammatic representation of the main process characterizing the attractive 
part of the three-nucleon interaction. The double line represents the excitation of a A 
resonsance (Ma = 1232 MeV) in the intermediate state. 

shaded region emphasizes the increasing importance of the three-nucleon contribution 
as the neutron drip line (^^O) is approached. Typically, in this case its size is such 
that {Vijk)/{vij) ^ 10%. For isospin symmetric nuclear matter, the role of the 
three-nucleon force in determinig the saturation properties is illustrated in the right 
panel of fig. (1.7). The density dependence predicted by several different equation of 
states are plotted both with (green and red lines) and without (black line) inclusion 
of three-nucleon interactions. 
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Figure 1.7. Left panel: Ground state energy of Oxygen isotopes. The red hne (ind) has 
been obtained using a chiral N^LO two- nucleon potential evolved through the similarity 
renormalization group technique. The results represented by the blue line (full) also 
include the contributions of a N^LO chiral three- nucleon potential. All calculations 
have been carried out within the SCGF formalism [?]. Close to the neutron drip line 
the effect of three- nucleon forces (the shaded light blue region labelled (5), turn out to 
be quite important. Right panel: Equation of state of isospin symmetric nuclear matter, 
obtained with and without inclusion of three-nucleon interactions [20]. 
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Chapter 2 

Correlated basis functions 
formalism 

In non-relativistic nuclear many-body theory (NMBT), a nuclear system is seen as a 
collection of point-like protons and neutrons whose dynamics are described by the 
hamiltonian: 

i j>i 

where t{i) and v{ij) denote the kinetic energy operator and the bare NN potential 
respectively, while the ellipses refers to the presence of additional many-body inter- 
actions. Due to the presence of the repulsive core, carrying out perturbation theory 
in the basis provided by the eigenstates of the non interacting system requires a 
renormalization of the NN potential. This is the foundation of many widely employed 
approaches developed to describe both finite nuclei and in infinite nuclear matter. 
For example, the schemes based on the no-core Shell Model [23], the Similarity 
Renormalization Group [24] and G-matrix perturbation theory [25], aim at obtain- 
ing the in medium NN scattering amplitude from a bare potential. Alternatively, 
the many-body Schroedinger equation can be solved using different approaches: 
self-consistent Green function theory [26], coupled cluster method [27], as well as 
stochastic or variational techniques. Typically, each one of these methods have a 
preferred environment, finite nuclei or infinite nuclear matter, where they have been 
applied successfully. However, stochastic and variational techniques have proven 
capable of providing accurate results both in light nuclei and uniform neutron and 
nuclear matter [28, 29, 30]. 

It has to be emphasized that within NMBT the interaction is completely de- 
termined by the analysis of exactly solvable two and three- nucleon systems. As 
a consequence, the uncertainties associated with any many-body calculations are 
decoupled from the the determination of the dynamical model. In principle, given 
the hamiltonian, the properties of nuclear systems ranging from deuteron to neutron 
stars can be obtained in a fully consistent fashion, without including any additional 
adjustable parameters parameters. 

In this work we use a scheme formally similar to standard perturbation theory, 
in which non perturbative effects arising from the short range NN repulsion are 
embodied in the basis functions, to be determined through a variational procedure. 
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2. Correlated basis functions formalism 



In the case of nuclear systems, the choice of the trial wave is particularly critical. To 
see this, consider, as an example, the fermion hard sphere system, i.e. a collection 
of particles interacting through a potential which is vanishing at all values of the 
interparticle distance r, except for the region r < a, in which it is infinite and 
positive. 

Using the eigenstate of the non-interacting Fermi gas, 
M Y(^j^^ ' "^FGi there is no way to handle the infinity, as the 

j behavior of the two-particle wave function is linear 

I ^FG near a nodal point (see picture on the left). As a 

consequence, a finite region in which V{x) and 
overlap is always present, and the expectation value 
{V) oo. On the other hand, a good trial wave 
function must reflect the features of the actual ground 
state, 

Hence, the two-particle wave function associated with must be vanishing at r < a, 
and different from zero at r > a. The better we choose ^ we, the closer the 
variational energy is to the eigenvalue of the hamiltonian. 




2.1 Correlated Basis Function Theory 

In the Correlated Basis Function (CBF) approach, the basis states are correlated waxe 
functions, obtained from the corresponding wave functions of the non interacting 
Fermi gas through the transformation [31, 32]: 

'""^ " (nFG|i^tF|nFG)i/2 ' 

where I^fg) is the determinant of single particle states describing N non-interacting 
nucleons. The operator F embodies the correlations among the particles induced by 
the NN potentaial. It is usually written in the form 

N 

F(1,...,7V)=5 n ^ ' 

j>i=i 

S being the symmetrization operator, needed to take into account the fact that, in 
general 

The structure of the two-body correlation function, fij, must reflect the complexity 
of NN potential. As stated in the previous chapter, the calculations discussed in 
this Thesis have been performed using the truncated Argonne vq potential. As a 
consequence, we chose the same operatorial structure for fij 

n=l 

As already pointed out, the model provides a fairly accurate description of the 
correlation structure of two-nucleon bound state. The shapes of the functions f^{rij) 
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are determined through the functional minimization of the expectation value of 
nuclear hamiltonian in the correlated ground state 

Eo = {0\H\0) . 

The evaluation of the above expectation value involves a degree of difficulty that 
rapidly increases with the number of particles in the system, A. It has been computed 
exactly using Monte Carlo techniques only for A < 16. For larger nuclei, as well 
as for infinite nuclear matter, approximate calculations are carried out using the 
cluster expansion formalism, to be discussed in the next Section. 



2.2 Cluster expansion formalism 

The correlation operator F defined above, must exhibit a factorization, or clustering^ 
property, dictated by the finite range of the NN interaction. Suppose that we pick 
a subset of p particles, labelled zi, . . . These particles are then moved away 
from the remaining ip+i, . . . ^in particles. Under these conditions, the A-particle 
correlation operator F must factorize into the product of two pieces, according to 

. . . , TV) ^ . . . , ip)FN-p{ip+q, . • • , ^tv) , 

implying that the two subsystems become completely independent of one another. 
The above property provides the basis of the cluster expansion formalism, that allows 
one to express any matrix element of a many-body operator between correlated states 
as a sum of contributions associated with isolated subsystems (clusters) consisting 
of an increasing number of particles. 

As an example, let us the consider the expectation value of the hamiltonian 
neglecting the three-nucleon potential. We will closely follow the derivation of the 
corresponding cluster expansion given in [32]. The starting point is the definition of 
the generating functional: 

I{/3) = (O|e^(^-To)|0) , (2.1) 
where |0) = . . . , A^)fg) is the actual correlated ground state of N particles and 

\Pi\<PF 

the kinetic energy operator. Starting from eq. (2.1) we can rewrite the expectation 
value of the hamiltonian in the form 



{H) = {0\H\0) = To+ ^ln/(/3) 



(2.2) 

/?=o 



Let us now expand the rhs of the above equation according to the cluster formaUsm. 
The functional in (2.1) involves the degrees of freedom of all the A'^ particles. As 
the particles are indistinguishable, any subset oi p < N particles can be arranged in 
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2. Correlated basis functions formalism 



N\/{N — p)\p\ ways. Defining a generating functional for each subset, we get: 
/i(/3) = (i|e/5Wi)--?)|i) , 



/i,...,n(/3) = m , 

where the indices label the non-interacting single-particle states of the 

"vacuum", or Fermi sea, is the corresponding kinetic energy, v{ij) = Vij and 

Let us now reexpress /(/3) in terms of n-body generating functionals (n < N)^ 
with = 1, i. We start noting that lij is close to the product of li and /j. 

It would in fact be exactly equal if we could neglect the interaction '^(12) and the 
associated correlation ^2(12). This observation suggests to rewrite lij as. 



IJ 5 



where the deviations of Yij from unity arise from correlation effects. Relabeling 
Yi = li we can describe the above functions in terms of the Y^^... according to 



implying 



W) 

lnI(/3) 



Yi , 

y.y.y. . 



(2.3) 



^ In + ^ In Fi,- + . . . + In Y^,,,^ . (2.4) 

i j>i 

Here the p-th term gathers all contributions involving all possible interaction among 
p particles. In a diagrammatic language, this term can be represented by a j)- vert ex 
diagram representing the nucleons in the cluster, connected by lines corresponding 
to dynamical and statical correlations. Substituting the above equation in eq. (2.2) 
we can write the ground state expectation value of the hamiltonian in the form 



(H) = To + {AE)i + {AE)2 + ... + {AE)n, 
where the contribution of the p-body cluster is given by the p-th term^, 

d 



(2.5) 



E 



ii<...<ip 



13=0 



In order to express eq.(2.5) in terms of the functions Ijj , we invert the relation in 
(2.3), 

Yi =Ii, 



Yij =(/i/j)-ili 



^As the particles are indistinguishable, the antisymmetrization of two body matrix elements only 
requires antisymmetrization of the state 
^Note that (A^)i = 0, since /, = = 1. 
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The resulting two-body cluster term is 



1^ 



dli a/j 



(2.6) 



As the number of particles in the cluster increases, the difficulties involved in 
the evaluation of the corresponding contribution become more and more severe, 
making numerical calculations impossible. However, the diagrams associated with 
the different cluster contributions can be classified according to their topological 
structure, and selected classes of diagrams can be summed up to all orders solving 
a set of integral equations, called Fermi Hyper-Netted Chain (FHNC) equation 
[33, 34]. 

2.3 Effective interaction 

The in medium effective interaction is usually defined starting from eq. (2.5) and 
setting 

(H) = {Ofg\To + FefflOre) • (2.7) 

The idea underlying the effective interaction approach is that the above equation 
yields screened potential, that is well behaved, and can therefore be used to carry out 
perturbative calculations in the basis of eigenstates of the non interacting system. 

Here, we will follow a procedure developed in Refs. [35, 36], whose authors 
derived the effective interaction from a cluster expansion of the left-hand side of eq. 
(2.7) truncated at the two-body level. 

Comparison between eqs. (2.7) and (2.6) shows that V^ff can be expressed in 
term of the effective two-nucleon interaction Wi2, defined through 



(OfgIK^sIOfg) 



Y.{^j\wi2\ij)c 

i<j 



(2. 



with 



/l2,[t(l)+t(2),/i2] + /l2t'(12)/l2 



2m 



/l2,[V2,/i2] + /l2t^(12)/i2 . 



(2.9) 



In the second line of the above equation we have exploited the fact that the correlation 
functions only depend on the relative distance between the interacting particles, r. 
Moreover, in the limit N oc, the normalization /ijl^^g of eq. (2.6) differs from 
unity terms 0{1/N) at most, and can therefore be neglected. The details of the 
calculations leading to the final expression of Wi2 are given in Appendix B. The 
result is 



1 



Wi2 



m 



(V/12) ■ V + t;(12)/i2) 



/l2(--(VVl2) 

/l2(-^(V2) + t,(12))/i2 = E<//WO^ 
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2. Correlated basis functions formalism 



Note that only the static part of the effective interaction is retained. The results of 
numerical calculations show that the contribution of the term (V/12) • V, yielding 
an expcitely momentum-dependent contribution through the exchange part of the 
matrix elements, is in fact very small compared the one arising from the static term. 

In order to include in V^ff the effects of three- nucleon forces, we follow the 
procedure adopted in Ref. [36], based on the approach originally proposed by 
Lagaris and Pandharipande [37]. Within this scheme, interactions involving three or 
more nucleons are taken into account through a density dependent modification of 
the NN potential at intermediate range, where two-pion exchange is believed to be 
the dominant mechanism. Neglecting, for simplicity, the charge-symmetry breaking 
components of the interaction, the resulting potential can be written in the form [37] 

vi^j)= E mrij) + vf{rij)e-^^^ + v-sirij)]0^^, (2.10) 

n=l,U 

where v'^^ and denote the long- (one-pion-exchange) , intermediate- and short- 
range part of the potential, respectively. The above modification results in a repulsive 
contribution to the binding energy of nuclear matter. The authors of Ref. [37] also 
include the additional attractive contribution 

AETNA = 72P'(3 - 2/3)e-^3P ^ (2.II) 

with /3 = {Pp — Pn) / {pp + pn) 1 where pp and pn denote the proton and neutron density, 
respectively. The values of the parameters 71, 72 and 73 appearing in Eqs.(2.10) and 
(2.11) have been determined in such a way as to reproduce the empirical binding 
energy and equilibrium density of nuclear matter [37]. 

Given the bare potential, the effective interaction is determined by the correlation 
operator /12. Its shape can be obtained from the functional minimization of energy 
at two-body cluster level, yielding a set of coupled differential equations, to be solved 
with the boundary conditions, 

'1 n= 1 



= \ n = 2,3,4 
fn{r>dt) = n = 5,6 



and 

= n = l,2,3,4 



dr 
dfn 
dr 



r=d 

= n = 5,6 

r=dt 



where d and dt are variational parameters. The above condition simply express the 
requirements that i) for relative distances larger than the interaction range, the 
two-nucleon wave function reduces to the one describing non interacting particles 
and ii) tensor interaction have a longer range, implying dt > d. We have solved the 
Euler-Lagrange equations for a wide range of nuclear matter density using the values 
for d and dt obtained from the highly refined variational calculation of Ref. [38], 
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carried out within the FHNC-SOC scheme. The derivation of the Euler Lagrange 
equations is discussed in Appendix C . 

In figure (2.1) the different components of the effective interaction at density 
p = 0.16 fm~^ are compared to the corresponding components of the Argonne vq 
potential. It clearly appears that the screening arising from correlations leads to a 
significant quenching of the bare NN interaction. 




Figure 2.1. Comparison between the components of the bare Argonne vq potential (blue 
dashed line) and those of the CBF effective potential (pink solid line) at p = 0.16 fm~^. 
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2. Correlated basis functions formalism 



2.4 Ground state energy 



Within the scheme described in the previous Sections, the final expression of the 
ground state energy reads 

E 

N ^^"^ ' 2 



where (Tq) is the Fermi gas kinetic energy and we have exphcitly written the direct 
and exchange contributions to the matrix elements of wi2. 

One of the attractive features of our approach is that it allows one to treat any 
forms of isospin asymmetric matter at the same level of accuracy. Labeling with 
the fraction of particles of type 7^, we can write 



7=1 



fP ^ 

7=1 



1 . 



Each set of uniquely defines a type of matter, for example, pure neutron matter 
corresponds to 

1 





Xi X2 



Xs X4 



The expression of the ground state energy of any kind of matter is (see appendix B) 



E 
N 



E 



2m 



7 ^ P 



x,x^ I d'r w^, [a;^ - B:;/{pF^,r)i{pF,^r)] (2.12 

7M n=l 



where is the component of the effective NN interaction associated with the 
operator O^, pf,j = (Qn'^Pj)^^^ and the Slater function £ is defined as 



while 



An 
^7M 



denote the matrix elements of the six operators in spin-isospin space, where. Their 
exphcit expressions for the direct channel are 



A 



with a{0) — (3cos^^ — 1), 6 being the angle between r and the z-axis. The 
corresponding expressions for the exchange channel are 

/ 1 1 1 1 

2 2 

2 2 

\ 4 





Oi 


O2 


O3 


O4 


O5 








/ 1 


1 


1 


1 


a{e) 


a{e) 




Pt pi 


1 


1 


-1 


-1 


-a{e) 


-a{9) 




p t ^ t 


1 


1 


1 


-1 


a{6) 


-a{e) 






I 1 


-1 


-1 


1 


-a{0) 


a{e) 


) 



B 



a{9) 


a{e) 


-a{6) 


-a{6) 





2a{0) 





-2a{e) 



^7 = 1, 2, 3, 4 corresponds to spin-up proton (p t), spin-down proton (p spin-up neutron (n t) 
and spin-down neutron (n respectively. 
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The matrices are obviously symmetric. Equation (2.12) is completely general, and 
allows one to calculate the ground state energy for any forms of nuclear matter, 
once the are fixed. In fig. 2.2, taken from Ref. [39], the energies of pure neutron 
matter (lower panel) and isospin symmetric nuclear matter (upper panel) obtained 
from eq. (2.12) are compared to the results of Refs. [40] and [41]. In [41] (solid 
lines) the calculation have been carried out using a variational approach based of 
the FHNC-SOC formalism, with a nuclear hamiltonian including the Argonne vi^ 
and Urbana IX potentials. The results of [40] (dashed line of the lower panel) is 
have been obtained within the Auxiliary Field Diffusion Monte Carlo (AFDMC) 
approach using, the truncated Argonne interaction and the same three-body 
potential. It clearly appears that our effective interaction provides a description in 
excellent agreement with that resulting from state-of-the-art theoretical approaches 
over a broad density range. Note that the correct saturation properties of symmetric 
nuclear matter are only obtained thanks to the inclusion of the effects of three- 
and many-body interactions. It should also be emphasized that, using eq. v^^ and 
the model of Ref. [37] , amounts to effectively including contributions of clusters 
involving more than two nucleons. 

As a final remark, it is worth mentioning that our approach, that does not 
involve any adjustable parameters, also provides a very reasonable estimate of the 
compressibility module^ of symmetric nuclear matter: K — 230 MeV. 



^ To be defined in eq. (4.9). 
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Figure 2.2. Upper panehEneigy per particle of symmetric matter as function of the 
density. Lower panel: Same as in the upper panel, but for pure neutron matter. The 
diamonds and the solid lines represent the result obtained from eq. (2.12) and those of 
Ref. [41], respsctively. The red dashed line in the lower panel shows the result of the 
AFDMC approach [40]. 
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Chapter 3 

Landau theory of normal Fermi 
liquids 

Landau theory provides a conceptual framework for understanding the main features 
of neutron matter near T ~ 0. The theory describes fermionic systems that are 
assumed to be in normal state, i.e. not showing any tendency to aggregate particles 
into clusters exhibiting bosonic features. In an infinite system, these clusters are the 
result of an attractive interactions, the occurrence of which completely distorts the 
physics near the ground state. The best known example is liquid ^He-*^, which is so far 
the only system in which quantum effects appear before solidification. Further exotic 
examples are electrons, both in the conductive band of metals and in white dwarfs, 
and, of course, nuclear matter. It is a remarkable fact that these systems, while being 
very different from one another, can actually be described by the same theory. This 
is due to the combined effects of the fermionic nature of the constituent particles and 
the zero temperature limit. In this context, statistical correlations largely dominate 
over dynamical correlations. The system is "frozen", almost degenerate, its state 
being very close to the corresponding state of the non-interacting system, where 
only statistical correlations induced by the Pauli principle matter. In the ground 
state, the particles fill up the lowest energy levels, up to the Fermi energy ep, while 
the states of higher energy are empty. As a consequence, the Fermi energy coincides 
with the chemical potential of the system, /x. 

Interacting systems close to the ground state, can be equally well described in term 
of two set of quantities, the single-particle energy levels, 65, and the corresponding 
occupation numbers, as described by the distribution function by n{es). The latter 
provides the the probability density for state of energy Cg to be occupied. In a cold, 
non interacting, Fermi gas this distribution function reduces to the step function 

n(e,) = ^(e, - ep)- ■ (3.1) 

In principle, switching on the interaction, the system can no longer be described in 
term of single particle states. In order to evaluate any macroscopic observables, one 
should find the stationary states of the system as a whole. Nevertheless, low-lying 
excited states are strongly influenced by the ground state itself, since interactions 
can only affect the states near the Fermi level, where empty and occupied states 

^Although the theory can also be extended to describe ^He-^He mixtures. 
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are very close to each other. Particles with energies lying deep in the Fermi sea still 
behave as free particles, being largely unaffected by interactions. As a consequence, 
all scattering properties, giving rise to multiparticle excitations, are restricted to a 
narrow region near the Fermi surface. 

The features of the energy spectrum of the system in this region of energy 
can be inferred from very general considerations, which are valid regardless of the 
magnitude and the specific nature of the interaction. Each low-lying excited state 
of a macroscopic system can be seen as simple combination of a set of elementary 
excitations, dubbed "quasiparticles", which largely behave like non-interacting 
particles. Their dispersion relation {uj — cj(k)), i.e. the equation linking energy 
end momentum, is the fundamental brick to build up the excitaton spectrum of the 
system. 

The concept of elementary excitations of a many-particle system can be best 
understood considering the vibrations of a lattice[42], discussed in the next Section. 

3.1 Quasiparticles 

Let us consider a perfect crystal, and focus on lattice excitations. As long as 
the amplitudes of its vibrations are small, the lattice be regarded as a set of 3A^ 
coupled harmonic oscillators {N denotes the number of atoms). Introducing normal 
coordinates, the system is reduced to a collection of 3A^ oscillators with frequencies 
uji^ and its energy takes the form 



where the ui are integer non-negative numbers. Lattice vibrations can be can be 
described as a superposition of monochromatic plane waves propagating through the 
crystal, characterized by a wave vector k and a frequency f^(k) which, in general, is 
not a single- valued function of k. For small k, we can have, for example, a linear 
dependence on the wave vector, i.e. fi(k) = u{9^ 

An alternate description is based on quantum mechanical wave-particle duality, 
stating that each plane wave can be associated with a particle carrying momentum 
k and energy ri(k). According to this picture, an excited state of the lattice can be 
seen as a collection of particles, called phonons, moving freely within the volume of 
the crystal, and its energy is determined by the number of phonons in the state i. 

In phonon language an anharmonic term in the energy can be interpreted as 
a phonon-phonon interaction, giving rise to scattering processes or creation of 
additional phonons. 

The description of lattice excitations in terms of phonons is closely related 
the description of normal Fermi liquids in terms of quasiparticles, proposed by 
Landau. At near zero temperatures, quasiparticles are the relevant degrees of freedom. 
However, in the presence of interactions (the anharmonic term, in phonon language), 
quasiparticles become unstable, as they do not describe the exact eigenstates of the 
system, and transitions via decay or scattering lead to a damping of the excitations. 
Decays occur for high energy excitations, while scattering becomes the dominant 
damping mechanism with growing number of excitations. Using simple phase-space 





3.1 Quasiparticles 
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considerations [43] , it can be shown that the quasiparticle concept is meaningful only 
in the vicinity of the Fermi surface and at low temperature. The inverse quasiparticle 
lifetime can be written as 

^ a(e - + , (3.3) 



T 



where fi is the chemical potential and a, h are positive constants. In the vicinity 
of the Fermi surface and a sufficiently low T, quasiparticles are long-lived states, 
and pretty much behave like stable particles. Equation (3.3) states the domain of 
applicability of Landau theory. Within this range all the ground state properties 
of the system can he described through the distribution function associated with the 
ground state of non interacting quasiparticles^ which in turn corresponds to the 
corresponding state of interacting physical particles. 

Adding an interaction, or raising the temperature, leads to a change in the 
distribution function near the Fermi surface. As quasiparticles are elementary 
excitations of the system as a whole, Landau theory assumes Sp to be a functional 
of quasiparticle distribution function, i.e. 

= '£^p(^p) • 

As a consequence, the total energy of the system can be written as, 

£; = £;o + ^5,(k)5n,(k) . 

k(7 

Note that, as in the phonon description of lattice excitations, the information on 
the ground state energy £"0 is totally lost. The definition of quasiparticle energy can 
be obtained performing a functional derivative of E 

"-^^^ - 9n.(k) • 

On the other hand, the distribution function ncr(k) at equilibrium is determined by 
minimizing the total energy 

5E{n) = T5S{n) + ^5N{n) . 

In the above equation, S{n) denotes the entropy, which is the same as in the free gas, 
its form being derived from purely combinatorial considerations^. As a consequence, 
one finds the ralations 

5£; = ^6,[n,(k)]5n,(k) , 57V = ^5n,(k) , 55 = - ^ 5n,(k)ln-^^^, 

ka ka ka ^ ^^^^^ 

leading to an implicit equation for the quasiparticle energy, as the Fermi-Dirac 
distribution 

n^(k) = {1 + e^{^M^)\-^}y^ ^ (3.4) 
in turn depends on £:cr[ncr(k)]. 



strong assumption that has to be done is that the number of interacting states is the same as 
the number of states in the non-interacting system. 
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3.2 Landau Parameters 

The quasiparticle energy is an additive quantity, as long as the number of quasi- 
particles is very small with respect to N. Under these conditions, their interaction 
probability is negligible and quasiparticles can be assumed to be free, the quasiparti- 
cle energy is defined as the first functional derivative with respect to Sn and terms 
0(5n^) can be neglected. 

Switching on quasiparticle interactions, or warming up the system, we need to 
also include the terms 0(5n^). As a result, we obtain 

SE = ^eO(k)5n,(k) + \ /aa'(k,k')<5n,(k)(5,Kk') , (3.5) 

kcr aa'kk' 

where is the energy of non interacting quasiparticle and /cro-K^^^O describes 
the interaction between two quasiparticles carrying momenta k and k'. Since 

(k) — ^s(k) — n^(k) is appreciably different from zero only near the Fermi 
surface, all the momenta involved in /crcr/(k,k^) are restricted to the Fermi sphere 
k ^ k' ^ kp. As a consequence, at fixed density, i.e. fixed kp, /o-o-K^os^) only 
depends now on the relative orientation of the two momenta. 

The function /crcr/(k,k') must reflect the main properties of the interaction 
between bare particles. In the case of electrons, for example, it must be spherically 
symmetric and involve both the direct and exchange terms. In the case of nuclear 
matter, on the other hand, the composite nature of nucleons results in the appearance 
of a non central component, non-diagonal in spin space, that is usually taken into 
account through a tensor term according to 

/,,,(cos^) = f{q)I + g{q){(7i • (72) + h{q)SM) , 

with q= (k-kO/2 and 

^12(qj = o -— ^ [CTi • (72) . 

The three functions f{q),g{q) and h{q) embody all features of NN interaction, 
within the operatorial representation employed in our work, i.e the Argonne v^. Note 
that here the tensor term in has its counterpart in momentum space. 

Landau parameters are generally defined as the projections of /^^/ on the basis of 
states of definite angular momentum. Expanding f^ra' (cos 6) in Legendre polynomials, 
and multiplying by the density of states at the Fermi surface, iV(0), one obtains the 
dimensionless quantities 

00 

F,,,(cos^) = N{O)f,,,{cos0) = ^F^^ P^(cos^). 

The set of parameters F[^' = F^^Gi^Hi are known as Landau parameters. They 
provide a link between the microscopic dynamics and macroscopic observables, and 
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require either phenomenological information or a local microscopic hamiltonian to 
be determined. While not being directly obtainable from Landau theory itself, the 
Landau parameters have very important implications. Depending on their values, 
the system exhibits different phases, such as giant resonances or pion condensation. 

In the case of a purely central interaction, the derivation of macroscopic observ- 
ables at equilibrium is well known [43]. As an example, let us consider a scalar probe 
inducing density fluctuations of the system. In coordinate space this amounts to 
compression and rarefaction of matter density, similar to the propagation of a sound 
wave. In momentum space, on the other hand, it corresponds to an oscillation of the 
Fermi surface, that "inhales and exhales". This radial perturbation can be written 
in the form S'kjr = Ak^?, independent of spin. In this case, the compressibility of the 
system is found to be 

Now, suppose to add an external magnetic field B, oriented along the z-axis, which 
couples with the particle spin through a term ^ (cr • B), where fi denotes the 
magnetic moment. A measure of the strength of the response is provided by the 
magnetic susceptibility, defined as 



75 



(3i 

B=0 



3 /] 

Ma/3 = ^nOoL^ + -MT I - 



where M is the total magnetization of the system and the indices a and /? are refer 
to the euclidean space. Within Landau theory x is determined by the different 
behavior of spin-up and spin-down particles, which oscillates with opposite of phase 
5\<ip — ±k]/ . Note, however, that in the presence of non central forces, the Fermi 
surface is no longer spherical. In addition, the quasiparticle effective charges may be 
modified, implying in turn that the effective magnetic moment may have components 
which are not scalar under rotation, i.e. it may have different values in different 
directions. The general form of /i, taking into account the effects of tensor forces is 
[44] 

3 ^ _ [ ^a^(3 ^ 

k2 3 

where fi is the usual magnetic moment, while fiT represents the magnitude of non- 
diagonal term. The experimental value of /j^t has not been determined yet. In Ref. 
[45], Arima and collaborators proposed the value |//r| ^ 0.049|/x|. 

The calculation of the susceptibility in the presence of non central interactions 
was first discussed by Haensel [46] , and then carried out in a very general fashion 
by the authors of Ref. [44], who used different forms of both the tensor effective 
interaction and the magnetic moment. These results are typically expanded in 
powers of jit- However, this quantity is still largely unknown, and seems to be 
very small. In a normal system its effect seems to be completely negligible, and 
has been ignored in our work. Within this approximation, the diagonal element of 
susceptibility reduces to 



1 + Go \ 8 1 + Go 1 + G2/5 
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Note that in the hmit Hi 0, we retrieve the usual expression, suitable for a purely 
central potential. 

Corrections to the mass of the particles arise from a different kind of oscillation. 
Consider a displacement of the Fermi sphere by a small momentum q, without 
changing its size and shape. The energy-density of the system is increased by 




Assuming that near the Fermi surface the quasiparticle energies can be written 
using an effective mass expansion, i.e. that £^(k.) ^ |kp/2m for |k| ^ fci?, one 
can also obtain SE from eq. (3.5). Equating the two results, and substituting 
7V(0) = kFTn^/n'^, leads to 




2m 



implying 
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Chapter 4 

Calculation of Landau 
parameters 



4.1 Quasiparticle interaction 



As pointed out in the previous Chapter, the quasiparticle interaction must embody 
all the relevant features of the bare interaction, which in this work is written using 
the six operators needed to take into account the spin-ispspin dependence and the 
presence of a non central component. Obviously, as our analysis is restricted to pure 
neutron matter, (r • r) = 1, and the number of operators can be reduced to three: 
/, (cTi • cr2) and Si2(q). As a consequence, one has to introduce three different sets 
of Landau parameters, associated with the direct, spin-exchange and tensor terms of 
the interaction, respectively. 

Within Landau theory, the contribution of second order in 6n to the energy shift 
6£ = £ — So can be written in the form 



2V 



J2 Yl ^ni{ki)5nj{kj 



(4.1) 



'^5 J k^jko 



X {ij\{[f{q)i] + [g{q){a^ ■ ^73)] + [h{q)SM)]}\ij) , (4.2) 



where the sum is extended to quasiparticle pairs in states ij living on the Fermi 
surface. 

In order to derive the Landau parameters from the effective interaction described 
in Chapter 2, we have to combine the right hand side of the above equation with 
the corresponding result obtained using V^ff and the Hartree-Fock approximation, 
eq. (2.12). Variation of the energy with respect to 5n = n — npG yields the second 
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order contribution 



5^CBF =^ Yl YSni{ki)6nj{kj 



'{^J\{ 



r3 



+ 



1 d3r(?/;3 + W4){1 - P<.e-^i-'-)(o-i • 0-2)] + 
" J d\{w^ + wq){1 - P^e-''^-') Si2{r)] } . (4.3) 



Note that, owing to to the presence of the spin-exchange operator Per = [1 — 
(cTi • cr2)]/4, the operators involved in the above equation, Pcrrl^ Pa{o-i • cr2), . . ., 
are different from those appearing in eq. (4.1). Therefore, it turns out to be more 
convenient considering the matrix elements, rather than the operators themselves. 
For example, in the case of a pair of particle with spins we require that 

5^£l S^Scbf 



where: 



f{q) in \i\ n) +9{q) in i^i • ^21 n) +h{q) in is^m n) 



2P2{cosO 



On the other hand, in the sector we find 



f{q) {U \I\ U) +9{q) in 1^1 • U) +Kq) {U |5i2(g)| tl), 



-2P2{cosC) 

where P2{cos^) is the second Legendre polynomial and cos^ is the angle between 
q and the z-axis. Using the above results one can construct the symmetric and 
antisymmetric branches according to the standard definition 

nq)^g{q) + 2h{q)P2{cosO , 

showing that inclusion of the tensor term results in a change of the interaction 
strength depending on the direction of the momentum transfer. This feature is 
illustrated in fig. (4.1), displaying the functions (left) and (right) obtained 
from the matrix elements of the effective interaction of Chapter 2. 

Carrying out the second functional derivatives of Scbf^i we can readily identify 
the Landau parameters from 

27r 



f{q) 
h{q) 



V 

27T 

V 

47r 



J d|r|r2 



2{wi + W2) - go{3ws + 3w4 + wi-\- W2) 



d\Y\ 2{w^ + W4) - gQ{-wi - W2 + W3 + W4] 



d|r|r^^2 



+ We 



^The last term in eq. (4.3) is the Fourier transform of the tensor operator 5*12 (q). A simple 
derivation is given in Ref. [47], while a more formal discussion can be found in Appendix D. 
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-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



qx ifm ^) qx ifm ^) 

Figure 4.1. Density plots of the functions (left) and (right) obtained from the 
matrix elements of the effective interaction of Chapter 2, as a function of the x and 
y components of the momentum transfer q. Note that the functions are azimuthal 
symmetric, and shown for Qz = plane. The non spherically symmetric dependence 
arising from the tensor term is clearly visible. 



with 

2 J-i Y 

4.2 Landau parameters from the CBF effective interac- 
tion 

The values of the dimensionless Landau parameters F^, Gi and obtained from 
the matrix elements of the CBF effective interaction are listed in Table 4.1 and 4.2, 
for ^ = 0, 1, 2. Their density-dependence is displayed in fig. 4.2. 



Go 0-08- 




P(fm ') p(fm ^) 



Figure 4.2. Left panel: density-dependence of the Landau parameters and Gi of pure 
neutron matter obtained from the matrix elements of the CBF effective interaction. 
Right panel: same as in the upper panel, ut for the Landau parameters H^. 

Within Landau theory, the parameters of Tables 4.1 and 4.2 provide a description 
of both the static and dynamical properties of pure neutron matter. In the literature. 
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p{fm ^) 


Fo 




F2 


Go 




G2 


U.U4 
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U.I 00 


U.ib9 


0.U94 


0.08 


-0.481 


-0.452 


-0.143 


0.844 


0.142 


0.138 


0.12 


-0.374 


-0.544 


-0.180 


0.881 


0.107 


0.150 


0.16 


-0.263 


-0.605 


-0.230 


0.914 


0.067 


0.155 


0.20 


-0.120 


-0.675 


-0.307 


0.978 


0.027 


0.150 


0.24 


0.004 


-0.726 


-0.351 


1.006 


-0.017 


0.136 


0.28 


0.092 


-0.734 


-0.390 


1.027 


-0.052 


0.121 


0.32 


0.184 


-0.763 


-0.412 


1.039 


-0.088 


0.102 



Table 4.1. Landau parameters and Gi of pure neutron matter obtained from the matrix 
elements of the CBF effective interaction (see left panel of fig. 4.2). 



p{fm ^) 


^^0 


Hi 


H2 


0.04 


0.060 


0.004 


-0.018 


0.08 


0.072 


0.019 


-0.016 


0.12 


0.056 


0.036 


0.002 


0.16 


0.046 


0.060 


0.025 


0.20 


0.038 


0.072 


0.041 


0.24 


0.028 


0.078 


0.052 


0.28 


0.025 


0.084 


0.060 


0.32 


0.019 


0.090 


0.069 



Table 4.2. Same as in Table 4.1, but for the Landau parameters (see right panel of fig. 
4.2). 



one can find limited sets of Landau parameters, obtained from phenomenology and 
semi-quantitative treatment [51, 52, 53]. Formally F^, G^, can be related to the in 
medium forward scattering amplitude. Nevertheless, as pointed out in the previous 
Chapter, their physical interpretation is deeply related to the static properties of 
matter. One of the major achievements of Landau theory is that once the Landau 
parameters has been fixed, all dynamical properties can be naturally derived from the 
static ones. 

In principle, the description of short and long range correlations should be carried 
out in a consistent fashion, using perturbation theory in the basis of correlated 
states. While static properties can be derived directly from the equation of state 
(EOS)^, the calculation of the dynamic responses involves non trivial difficulties. A 
prominent issue is the treatment of collective modes, whose contribution becomes 
more and more important as the momentum transfer decreases. In particular, the 
calculation of the response function at |q| ^0 within the CBF formalism requires 
the description of the propagation of correlated one particle-one hole states through 
the medium. A consistent description of the response of isospin symmetric nuclear 



^In general, the EOS is a non trivial relation linking the thermodynamic variables specifying the 
state of a physical system. 
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matter, performed using the CBF effective interaction, has been reported in ref. [48]. 
However the analysis is restricted to the density channel. 

In this Thesis, we adopt a different approach. The dynamic response is no longer 
obtained microscopically, summing the contribution of all particle- hole pairs. We 
take instead the point of view Landau theory, in which the degrees of freedom are the 
quasiparticles, representing the elementary excitations of the system. This approach 
is both simpler and more general than that discussed in Ref. [48] , as it is capable to 
treat the coherent and incoherent contributions to the response on equal footing. 
All the main features of CBF interaction, reflecting both the nature of the bare 
interaction and correlation effects, are summarized by few parameters, the values of 
which must satisfy a set of stability conditions. 

The main effect taken into account by the CBF effective interaction is the 
screening of the nuclear medium, arising from short range correlations, that leads to 
a strong suppression of the repulsive core of the bare NN potential. This pattern is 
exactly what we observe in the Landau parameters. The large value of Go, compared 
to Fq, reflects the negative value of f^^ < 0. This means that the effective interaction 
between particles with opposite spin is attractive, thanks to Pauli principle that 
keeps particles with parallel spin apart. In a different system, such as ^He, this 
lead to the opposite relation, Fq ^ Go, since the interaction is not screened, and 
particles with opposite spin experience the repulsive short-range interaction. As a 
consequence /^^ becomes more repulsive than /^^. Within the CBF scheme, the 
repulsive core is screened, and the attractive interaction emerges. Also note that at 
low density the channel is the most important, and greatly influence the values 
of Landau parameters. 

Increasing the density (or the pressure) Fq increases more rapidly than the other 
parameters, as it is strongly related to the compressibility of the system. Its growth 
is not as fast as in ^He, where the rapid variation is again attributed to the repulsive 
core of the potential » \f^^\). At p ^ 0.16fm~^, the equilibrium density of 

symmetric nuclear matter, fig. 4.2 shows that Gq ~ 0, or \f^^\ ^ implying that 
the repulsive interaction between particles with parallel spin has the same magnitude 
as the attractive interaction. Increasing the density |/^^| becomes larger, in spite of 
screening. 

Let us consider the harmonic ^ = 1. The values of the Landau parameters are 
almost constant in the range shown in Fig. 4.2, since the centrifugal barrier prevents 
the particles from coming close to one another. In particular, Fi < has a direct 
physical interpretation, since it is related to the effective mass. 

For incoherent excitations the stability conditions F^, G^ > —(2^+1) are satisfied. 
Nevertheless, in the low density limit the attraction in the channel appears to 
be dominant, and may lead to the appearance of an instability region, in which 
Fo < -1 for p < 0.04 fm-^ 

We now turn to the discussion of collective modes. The occurrence of these 
excitations depends on several relations between Landau parameters that were 
extended to the case of nuclear matter, with the inclusion of non central interactions 
at the end of 70s [54]. Comparing the time scale of collective modes with the 
quasiparticlep mean free path, (A), we can distinguish al least two different classes. 
When the length scale is longer than A, hydrodynamics may be used. In this limit 
the only relevant variables are density fluctuations and the average fluid velocity. 
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that correspond to ^ = and i = 1 distortions of the Fermi surface. The distortions 
arising from all higher harmonics are removed by collisions. In the opposite limit, 
i.e. short length scale, the higher harmonics persist, because of the absence of 
collisions, and a detailed description in term of quasiparticle distribution function is 
required [42]. The range — 1 < Fq < corresponds to strong Landau damping in 
the colhsionless hmit, since vc < vf-, where vc denotes the velocity of the collective 
mode. The collective mode are destroyed by incoherent excitation. Nevertheless, we 
have undamped oscillation in the hydrodynamic limit. For small positive value of 
Fo > the occurrence of the collective mode depends on the magnitude of Fq. A 
sufficient condition for stability was first found by Haensel, who studied the positive 
solutions of the the zero sound dispersion relation [55]. Considering the harmonics 
^ = 0, 1 he obtained 

Another sufficient condition was derived by Mermin [56] : 

y > , for£ = 0,l , -^V>— T^- (4-5) 

However, the above relation is more restrictive than the forst one for Fo > 0^. Note 
that with F we denote either F or G. In the same paper [56], Mermin demonstrates 
one more relation between Landau parameters^: 

I 2i+l ^ ^ 2i+l / 

In general, eq. (4.6) implies that at least one set parameters F, G satisfies condition 
(4.5). Hence, at T = at least one of the two channels must exhibit a zero sound 
wave. The values reported in fig. (4.2) satisfy the condition (4.5): zero sound always 
occur in the spin channel. On the other hand, in the density channel the condition 
of eq. (4.4) is fulfilled for p > 0.32 only. 



4.3 Static properties of neutron matter 

Before discussing the application of Landau theory to the calculation of the dynamic 
structure functions, in this Section we report the results of the calculations of a 
variety of equilibrium properties of neutron matter at T = 0. 

We will focus on effective mass, m*, isothermal compressibility, x^, ^.nd magnetic 
susceptibility, x^i which, in the static limit, can be related to the energy-density, 
density and spin-density responses, respectively. As these quantities can be obtained 
both from matrix elements of the CBF effective interaction and from the Landau 
parameters listed in Tables 4.1 and 4.2, the analysis discussed in this Section provides 
a valuable consistency test of our approach. 

^We remind the reader that for Fq < solutions can exist, but are damped, as < vf 
^ While this is true for neutron matter, for nuclear matter there are different isospin states and a 
general condition involves all channels [57, 58]. 
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Single-Particle Spectrum and effective mass 

Within the CBF approach, starting from the expression of the energy 



p,cr 



^ = E £,Mp) + ^ E E [d'v\A{r) - B{r)e 



p,p' aa' 



^ip r -ip r 



na{p)n,'{p') . (4.7) 



Ej>i(U>12^(r)|u)a 

we can define, the single-particle spectrum in Hartree-Fock approximation as 



2m 



j>i 



P 

2m 



p' cr 

The effective mass is then defined as 

m* = p 



+ / ^EE / d'r[A(r)-B(r)e^P-e-^P^-]n.,(pO • (4- 



dea{p) 



dp 



Within the CBF effective interaction approach, the potential is perfectly known 
and well behaved, and we can apply the above expression without any problems. 
Hence, the effective mass is defined for all values of the momentum, p^. On the other 
hand, starting from eq. (4.7) we can vary the distribution function to obtain an 
estimate of the energy in the Landau scheme. Note that this is a general procedure, 
in which the interaction is unknown, and can be different in different systems. The 
most convenient procedure consists in changing variable, so that the derivative is 
applied to 5n, which is different from zero only near p = p^^ The inset of fig (4.3) 



^ The derivative with respsct to p is performed as follows 
Vp(7(p) = Vp^(y(0)-y(p-p'))n(p') 

= -E^P^(P - P')"(P') = - j d^'vBir) (Vpe'P-) ^e'P'-n(p') 

P' P' 

^In this case, we perform a change of variable and then use the Landau ansatz dn/dep 

q=p-p' 



Vp[/(p) = Vp^(y(0)-y(p-p'))n(p')= <^ = 

P' 

= -^l/(q)Vpn(p-q) = - / d^rBW^nq)-^"^^"''^ 

q q 



VpCp-q 
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the ratio between effective and bare neutron mass, {rrf /m)^ evaluated p — pp^ 
as a function of matter density. It clearly appears that the results obtained from 
the CBF single particle energies (dots) and from Landau theory (solid line) are the 
same within the numerical accuracy. 



Compressibility 

Compressibility provides a measure of the "elasticity" of a system, determining its 
response to a density compression. Classically, it is related to the speed of sound 
Vg ^ p at equilibrium. From the definition 

through standard thermodynamic relations one obtains 



1 



— + 2 — 
^ dp^ dp 



£ 



where £ — E/N denotes the energy per particle^. 

Using the above definition we can easily evaluate K from eq. (4.7) and compare 
it with the result obtained from the corresponding definition of Landau theory, eq. 
(3.7). The comparison shown in the left panel of fig. (4.3), indicates that the results 
of the two approaches are in perfect agreement with one another over a broad density 
range. 



Spin susceptibility 

Another important thermodynamical property is the spin susceptibility, that mea- 
sures the response to an external magnetic field under the assumption that orbital 
effects can be safely ignored (see eq. 3.8). Here, we outhne the standard procedure to 
evaluate the spin susceptibility within the CBF framework. The energy per particle 
of asymmetric nuclear matter can be accurately approximated using the expression 

^ E{a, /3, 7) = ^0 + E^a^ + ^^/^^ + E^rl^ , (4.10) 

with 

a ^ {x'i - X4) + (xi - X2) 

P = (X3 + X4) - (xi + X2) (4.11) 
7 = (X3 - X4) - (xi - X2) . 



^In studies of isospin symmetric nuclear matter it is useful to introduce the so called "compress- 
ibility modulus", defined as 



(4.9) 



From the above equation, it follows that Koo ~ at saturation density (po = 0.16 fm~^), where the 
first derivative of 8 vanishes. Note that interpretation of Kog is strictly associated with symmetric 
matter, since the energy of pure neutron matter does not exhibit a minimum. 
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Figure 4.3. Left panel: Compressibility of neutron matter, normalized to its Fermi gas 
value, as function of density in units of po = 0.16 fm~^. The dots and the solid line 
correspond to the results obtained from Landau's theory, eq.(3.7), and the equation of 
state computed using the CBF effective interaction, eq.(4.7), respectively. Right panel: 
same as in the left panel, but for for the spin susceptibility, - The dashed line shows 
the susceptibility obtained from the dynamic spin structure function including only the 
incoherent contribution. The inset of the left panel shows the density dependence of the 
ratio between effective and bare neutron mass. 



The above coefficients define the spin-isospin content of matter: for example, x\ — l/4: 
for all values of A yields E/N — £^snm = while pure neutron matter corresponds 
to xi = X2 = and X3 = X4 = 1/2, E/N £^pnm = + Er- Note that this 
obviously implies that Ej- can be identified with the symmetry energy. 

Let us consider fully spin-polarized neutron matter. The two degenerate ground 
states, corresponding to X3 = 1 and X4 = (o: = 1, spin-up) and X3 = and x/^ — l 
(a = — 1, spin-down), respectively, have energy 



E^ = E^ ^ £^PNM + E(j 



(4.12) 



with = Efj + EfjT-. For arbitrary polarization, a, the energy can be obtained from 
the expansion 



E{a) = E{0) + 



dE 
da 



1 d^E 



+ ... . 



(4.13) 



a=0 



As E must be an even function of a (see eq.(4.12)), the linear term in the above 
expansion must vanish. Hence, neglecting terms of order a^, we can write 



AE = E{a) - E{0) 



1 d^E 



a 



(4.14) 



a=0 



2 da^ 

In the presence of a uniform magnetic field B the energy of the system becomes 

Esia) = E{a) - ajuB, (4.15) 



where B denotes the magnitude of the external field, the direction of which is chosen 
along the spin quantization axis, and fi is the neutron magnetic moment. 



42 



4. Calculation of Landau parameters 



Assuming that equilibrium corresponds to a = ao^ i.e. that 

-/iS = , (4.16) 



dE 
da 



a=ao 

we obtain 

(4.17) 

From the definitions of the total magnetization 

M = /i(p3 - pa) = paop = /x^ I I , (4.18) 




a=0 



and of the spin susceptibility x 
we finally obtain 



M = xB , (4.19) 



In the right panel of fig. 4.3 we compare the susceptibility calculated with eqs. 
(3.8) and (4.20), normalized to the Fermi gas result. As for the compressibility, the 
two approaches are in very close agreement. It has to be pointed out that in the 
static limit the contribution of the Landau parameters associated with the tensor 
interaction, (listed in Table 4.2) turns out to be less than 0.1%, and their effects 
are hardly visible. Therefore, we can further simplify the calculations, neglecting 
non-central contributions in spin channel altogether. This approximation does not 
spoil the level of accuracy of our result. The dashed line in fig. (4.3) shows the 
contribution of single quasiparticle excitation (see section 5.4 below). This means 
that in the static limit a large fraction of the response comes from the excitation of 
a collective mode. 

In the left panel of Fig. 4.4, the magnetic susceptibility of neutron matter 
computed within Landau theory is compared to the results of Ref . [40] , obtained 
using the Auxiliary Field Diffusion Monte Carlo approach and nuclear hamiltonians 
including the truncated Vq and forms of the Argonne vis potential, supplemented 
with the Urbana IX three-nucleon potential. In the right panel, we compare the 
density dependence of the compressibility resulting from our calculations to the 
results of a variational calculations , carried out using the full Argonne vis NN 
potential and the Urbana IX three-nucleon potential. Our results appear to be in 
reasonable agreement with those obtained from highly refined theoretical approaches, 
the differences at large density being likely to be ascribable to the different treatment 
of three-nucleon forces, which are known to play a critical role at p > po- 
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Figure 4.4. Left panel: comparison between the compressibility computed within Landau's 
theory and the corresponding results obtained using the variational FHNC- SOC approach. 
Right panel: comparison between the spin susceptibility computed within Landau's 
theory and the corresponding results obtained using the Auxiliary Field Diffusion Monte 
Carlo approach. 
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Chapter 5 



Weak response of neutron 
matter 



This Chapter is devoted to the analysis of the Boltzmann-Landau (BL) equation, 
the solution of which can be interpreted as the dynamical form factor within Landau 
theory. Our discussion will closely follow the treatment developed in classic textbooks 
of quantum many-body theory [59, 42]. 

5.1 Response in the low-momentum transfer regime 

In the second quantization formalism, linear-response theory [49] provides a micro- 
scopic description of the dynamic form factor of interacting many-body systems. 
Consider an external scalar probe described by the potential 



with 77 = 0^, coupled to the density of the system, assumed to be in its ground state. 



The induced oscillation in the local density is proportional to the initial perturbation 
according to 



In the above equation, pq is the spatial Fourier transform of particle density operator, 
(Pq)nO denotes its matrix element between the ground state, |0), and an excited 
state and cjno — — E^) is the excitation energy. The operator p\ acting on 
the ground state generates one or more particle-hole (ph) excitations, or a collective 
mode (zero sound). 

It can be shown that in the limit of low-momentum transfer, q ^ 0, the matrix 
elements (pq)nO involving a single ph excitation do not depend on |q|, while the 



t/(r,t) =e^[^-^-(^+^^Wt/(q,cj) 



via 




5p(q,cj) = K(q,cj)t/(q,cj), 
where K is the response function, that can be written in the form 
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contributions of the collective mode and multipair excitations exhibit a oc V^q[| 
and oc |q| behavior, respectively. As a consequence, the |q|-dependence of the 
corresponding contributions to the response is(x|q|^, oc|q|^ and oc |qp. 

The formalism based on the BL equation does not take into account multipair 
excitation, but allows one to treat coherent and incoherent ph excitations on equal 
footing. Therefore, it turns out to be useful writing the response in the form 

XLandau 



E'lK)n0| 

n/0 



2WnO 



~^Xmultipair 



nO 



where the sum includes ph excitations only. In view of its |q|-dependence, neglecting 
Xmultipair IS cxpcctcd to be a reasonable approximation in the |q| ^ limit. 

As |q| 0, the main contribution associated with a single ph excitation exhibits 
a specific form, that can be easily understood. Consider the Feynman diagrams of fig. 
5.1, showing all possible processes involving a particle (t) with momentum pi and a 
hole (D with momentum up to second order in the interaction, represented by the 
dashed line. The second order process can be can be classified in three topologically 
different sets, usually labelled according to the Mandelstam variables 5, u, t, which 
reflect the different momenta associated with the loop. The fermionic ph loop, that 



order 



>•■•■< 



2nd Q^^g^ 



P2-q 



^ 



Pi 



7 



> 

Pi 




P2 



pi + q 




Figure 5.1. Feynman diagrams describing one particle-one hole interactions up to second 
order in perturbation theory. The second order diagrams in the second line are labelled 
by the Mandelstam variables s, t. In the limit |pi — P2I = ^ ^ only the loop in 
the last diagram is divergent [42, 60]. 



can be readily evaluated for the non interacting system, is referred to as Lindhart 
function. It turns out to be exactly the same as the function iloo obtained from 
the kinetic equation (see below). In the low momentum transfer limit, the first two 
second order diagrams can be disregarded with respect to the third, t-channel, one. 
In this last term, the integral of the fermionic loop is in fact divergent as g ^ 0. 
Summing up the t-channel diagrams to all order, one finds the expression of the 
response in the so called ring approximation 

ringappr ^ ^00 W -.x 

The same expression, with an interaction potential V expressed in terms of Landau 
paraeters is obtained from the solution of BL equation. 



5.2 Response near equilibrium 
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5.2 Response near equilibrium 

One of the main features of Landau theory is that, under certain assumptions, the 
dynamic form factors of the system are fully constrained by the Landau parameters. 
This is true if the excited state of the system, produced by the interaction with the ex- 
ternal probe, is such that the quasiparticle distribution in not far from that describing 
the system in equilibrium. In this case the non equilibrium, non homogeneous state 
can be specified by the quasiparticle distribution nka{'^,t) depending on position 
and time. Note that we can fix k and r at the same time, since the characteristic 
length of the spatial inhomogeneity. A, is much larger than the localization length of 
the quasiparticle A/. Typically the change in the distribution function occurs on the 
Fermi surface and is proportional to the temperature, implying Ap T/vf and due 
to the uncertainty principle, \f ^ Vf /T . In the region A » A/ the system is locally 
homogenous (globally quasi homogenous) and we can use a classic description based 
on a local distribution function to specify the momentum density at position r and 
time t. 

The expressions of all observables of the system must also exhibit an explicit r 
and t dependence. For example. 



where E{r^t) is the total energy density and /crcr^k, k^, r, r') describes the effective 
interactions between quasiparticles at a relative distance r — r^ At zero-th order 
the system is homogenous, implying 6^^(r,t) = ^kcr- "^^^ effective interaction is 
different from zero only if the quasiparticles are very close, typically at distances 
comparable with the scattering length, i.e. |r — ^ a ^ ^/^f- Since rikvK^'^) 
almost constant over a distance A » a, we can expand 5nkv^ t) — ^^kV^ (r, + . . . 
and disregard terms 0{a/X). It follows that 



5.3 Kinetic Equation 

Within Landau theory [43], all the relevant macroscopic quantities are linked to 
the miscroscopic dynamics through the quasiparticle distribution function. The 
evolution of the system is then determined by the kinetic equation for 6n. Let us 
start from the Boltzmann equation for quasiparticle distribution function 




6E{r,t) 



1 ^11// 



kcr kk Vcr^ 



with 




^^k,c7(r,t) 
dt 




(5.2) 
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where d/dt is a total time derivative and the quantity /[n], cahed cohision integral, 
describes the fluctuation of the number of quasiparticles in a phase space volume 
arising from their mutual interactions. The quasiparticle energy density is treated as 
the classic hamiltonian of the single quasiparticle, and the external potential ?7(r, t) 
is added to account for the interaction between the quasiparticle and the probe. 
Eq (5.2) can be expanded according to 



anka(r, ^) 1 f 5^ka(r, t) de\^^ dni^^{v, t) d{e\^^{v, t) + U (r, t)) 



} = I[n] 



dt h { 5kc^ dka 

with 

_ d{e + U) 

dr ' 5k ' 

and linearized, since the total distribution function is well defined only for |k| ~ kp. 
For small deviations from equilibrium 

5nka(r, t) = nka(r, t) - n^^ , 

where the n^^ is the distribution at equilibrium, is in fact non vanishing only in the 
vicinity of the Fermi surface. Therefore, the final expression is 

d6ni^a{r,t) ^ d6ni^a{r,t) de^^ _ dn^^jr, t) a(gka(r, t) + U (r, t) ^ 

dt dVa d\<ia d\<ia 

with 

dk^ dei^ dk^ dk ^ ' ' 



leading to 
dt 



+ v/e • VrSn\^a{r,t) 



+ |^v,- (^F(r,t)-gy,,Kk,kOVr5nkvK = , (5-5) 



where F{r,t) = -Vrf7(r,t) and 



_^ = -SisF - 6°^) . (5.6) 



5.3.1 The collision integral 

The evaluation of the collision integral is essential for deriving the dynamic response. 
The behavior of the system depends mainly on the value of ur where u is the 
frequence of the external perturbation and r is the quasiparticle collision time. 
According to Landau theory, the probability of collision between two excitations 
(quasiparticles) decreases according to the square of the spread of Fermi distribution 
function, implying in turn hence r ^ T~^. At large T, cjt << 1 and the perturbation 
propagates according to ordinary hydrodynamics, since locally the system is in 
thermodynamic equilibrium thanks to the many quasiparticle collisions. For example, 
the speed of sound is related to compressibility through the usual thermodynamic 
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relation u — ^J^P|^p^ and its damping depends linearly ovlt ^ [42] . On the other 
hand, if cjt ~ 1 diffusion is extremely damped, and the perturbation does not 
propagate. 

In the region of T ~ 0, i.e. cjt » 1, eq. (5.5) turns out to allow a new solution. 
Compared to the classical sound wave, this solution, driven by quantum effects, 
describes a somewhat different physics: a collective mode propagating in the liquid 
without experiencing any damping. This mode was first predicted by Landau, who 
dubbed it "zero sound". Let us assume that 

8n - e-^^\ (5.7) 

where l/cj is the time scale of local dynamic, which implies ijJt ^ \. It follows that 
the integral /[n], being of order of ^n/r, can be neglected compared to dn/dt. 

5.3.2 Spin-dependent solution in collisionless limit 

Consider an external magnetic field coupled to the spin according to U{r,t) = 
—gfiB^ • B(r,t), where B is aligned in the z-direction. The linear spin-density 
response is defined trough the magnetization of the matter 



dB 



/5 



(5.8) 



B=0 



where S = |B|, the magnetization M is given by 



M. = -^(% - %) = Ei^n^t - Sn^l) = ^ E '^^ ' (5-9) 

k k 



and the spin susceptibility is^ 

x^f:f(EK)i. (5.11) 

The kinetic equation will help us to determine the change in the Fermi distribution 
function (5n" in the direction of momentum k due to the external field H(r,t). 
Performing a Fourier Transform of eq. (5.5) we find 



- Vfc • q)(5nka(w,q) 



__ _ __ It; J 

kV^ 



. (5.12) 



^The attenuation coefficient is 7 ^ uj^t]/ pu^ where 77 is the shear viscosity, p is the density and 
u ^ Vf the quasiparticle velocity. Note that rj/p ^ v^r where Vf does not depend on T 
^In the hterature one can also find the slightly different definition 
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This above equation is similar to an integral equation, and can be solved by iteration. 
Substitution of the zero-th order expression in the last term yields 



v/c • q 



{w -Wk' q) del. 
v/c • q 



kt 



(w - Vfc • q) 



2 

5MB 



B + 



B + 



(5.13) 



implying ^n^^ = ~^^k^- It follows that in each point of the Fermi surface, corre- 
sponding to a momentum k, the oscillation is ^ (v/^ • q)/(it; — v/^ • q). 
The equation for 5n^ is ^ 



+ 2(q-Vfc) ^ 



S(q,^) + ^r(k,k05<,(q, 
k' 



w) 



We now make fro the solution the ansatz 



(5.14) 
(5.15) 

(5.16) 



Note that the above expression is in general quite reasonable: the first factor 
constrains the change in the distribution function to the Fermi surface, while z/^ 
provides a quantitative description of this variation in the direction of k. In the case 
of a 0-function distribution, one can rewrite the above expression using the angles 
6\ (j)' between q and k^ 

5nk/ = Q{kf + z/k/(0', (j)') - k') - Q{kf - k') ^ hvfjy^^{e\ 00^(4^ - £f) . (5.17) 

Here, we only consider azimuthally symmetric solutions, corresponding to z^kK^)^- 
Expanding both u\^f{cos9') and {cos 9, cos 9') in Legendre polynomials 

z/k^ = ^Pi{cos9')ui /^(cos6>,cos6>0 ^^Pi{cos9')f^ {cos 9), (5.18) 
i i 



eq (5.14) becomes 
cos6> 



cos9 



Em/ ^N^^k^ / n/\ o ^^^^ 



COS 9 



B{q,uj)y (5.19) 



We have used the the relations 



fn^n^ + fn^rii = (/tt - = {Sn^ - Sni) = /''(k,k )^nk/ 



-r(k,kO^<.. 



^In principle, eq. (5.14) can describe any kind of collective modes, that differ in the angular 
dependence vi^iO^cj)) as well as in velocity. The arbitrary dependence of v-^' on (0' ^(j)') determines 
different solution. For example, a e™"^ dependence corresponds to collective modes leaving the total 
volume of the Fermi sphere unchanged. 
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where the sum yields 



r r/^k' 

2^r(cosO^^'(cos^') = -2V J _r(cose)<5(£/ -£0,)^k'(cos( 



-V 



N(0) 
Air 



d(cos e')d(f)'f {cos Oi^k(cos 9') 



We can simplify the expression using the addition theorem of Legendre polynomials 
and expanding again in Legendre polynomials in space: 



2^ + 1 



^^k' I 

where 5 = (jj/{\q_\vf) and = VN{0)f^. The equation in (5.19) become: 



2£ + 



with 



J Zi jj, s 



The explicit expressions for ^ = 0, 1, 2 are the following: 



(5.20) 



(5.21) 



(5.22) 



flnn = 1 + - In 



s + 1 



^20^ 7: + 



3s' 



00 



3s' 



20 



(5.23) 



Considering only terms with £ = 0, 1, 2 in (5.21) we obtain three coupled linear 
equations for i/q, i^i and f2, that can be cast in the matrix form 



^ 1 + f^ooGo ^floiGi 

^iqGo ^(l + f2iiGi) ^f2i2G2 
y 020^0 3021^1 5(1 + ^^22^2) 

Finally, the spin-density response defined by 

1 



^ f 

1^1 I = -{9IJ'bB{ci,uj)) 
\ 1^2 



dnj 
del 



reads 



Sp"" = yY.^{ei - £/)N + Pi(cos^)i/i + P2(cos^)i/2] 



2V 



\ f^20 



(5.24) 



(5.25) 



Note that for £ ^ the integral of Legendre polynomials vanishes, since they are 
orthogonal. The value of vq can be obtained from the matrix equation above. Using 
only £ — 0,1 the matrix reduces to dimension 2x2, and 



X 



N{0) 



(5.26) 
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5.3.3 Coherent and incoherent contributions 

One of the striking properties of Landau theory is that it is able to manage coherent 
and incoherent response on equal footing. In eq (5.19) we can distinguish at least 
two cases depending on the value of the parameter s. 

s < 1 The speed of the perturbation, Ug = ^/|q|, is smaller then Vf^. Note that 
for CL;/|q| < Vf the denominator of eq.(5.19) can vanish, thus indicating the 
presence of a resonance between the external perturbation and the single quasi 
particle excitation. In this case the perturbation can lead to the excitation 
of incoherent single ph pairs, and is consequently strongly damped. This 
phenomenon is reminiscent of the propagation of light in a medium, in which 
the imaginary part of the refraction index is proportional to the magnitude of 
the damping. In the case under consideration, the same role is played by the 
imaginary part of the response, that can be obtained from 



00 



2 s-1 



1 - - In 
2 



s + l 



s < 1 



(5.27) 



Typically, at T = an incoherent qp excitation from the ground state consists 
of the disappearance of a particle with |k| < and the appearance of one 
with |k + q| > kp. The corresponding excitation energy is 



AE — Sk+q 



2m* 



1 



2m* 



2m= 



■(q' + 2k-q) 



(5.28) 



For foxed |k| and |q|, the maximum and minimum values of the above energy 
are AEmax = (q^ + 2|k||q|)/2m* and AEmin = (q^ - 2|k||q|)/2m*. These 
relations correspond to two hyperboles, delimiting the spectrum of incoherent 
excitation (fig. 5.2). For 5 < 1 the dispersion relation of external perturbation 
corresponds to the region between the two hyperboles, and the attenuation 
arising from the interaction with single ph excitations would made it vanish. 



Tm x{q ^ 0,w) 

incoherent ebccitations 



Fermi Gas 
>1 
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Figure 5.2. Response in the |q| ^ limit, as a function of the energy transfer uj. Incoherent, 
coherent and multipairs contributions are clearly visible [59] . 



5 > 1 The velocity of the wave associated with the external perturbation is larger 
than the Fermi velocity. Hence, it can only couple to the collective mode 

^Note that the condition cos^s = Us/vf in eq. (5.19) is the same as the condition for emission of 
Cherenkov radiation. However, here it refers to zero sound waves, triggered by a single quasiparticle 
excitation. 
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only. Generally, in this case the response is obtained adding an infinitesimal 
imaginary part to the frequency, i.e. replacing cj ^ cj — irf. From the formal 
point of view, this is the rigorous way to perform this calculation, letting 
77 ^ at the end. This replacement is irrelevant when dealing with incoherent 
excitations, since in this case the response has a non vanishing imaginary part, 
as shown by eq. (5.27) above. However, it turns out to be important in the 
region corresponding to 5 > 1, where the response can be written in the form 

= V[...] + iTT ^noois) S{1 + [Go + Y^^Poois)). 

A 

To obtain the above equaiotn, we have used the well-known expression to 
distinguish the real part, associated with the principal part integration, and 
the imaginary part yielding a Dirac 5-function. Let us now focus on the 
imaginary part alone. After some manipulations we find 

^oo{so)S{l + [Go + s^A]noo{s)) (5.29) 

6{s - So) 



-1 



\2soAnoo{so) + (Go + 5gA)f^[)o(5o 



where sq is the root if the equation obtained requiring the argument of the 
5-function to vanish. 

In this case the energy of the coherent wave is above the incoherent excitation 
spectrum. Therefore, the collective mode is undamped. We note that this 
excitation is propagating in Fermi liquids at T = 0. 

In fig. (5.3) we compare the imaginary part of the response function in both 
the density and spin channels. According to standard linear response theory the 
correlation function is proportional to the imaginary part of the response 

(cr(cj, -q)cr(0,q)) - Imx . 

The calculation has been carried out using the set of Landau parameters listed in 
the tables 4.1 and 4.2, including contributions associated with angular momentum 
components up to L = 2. For comparison, we also show the results of Ref. [61] 
(dashed line). It clearly appears that including the ^ = 2 parameters hardly affects 
the response in both the spin and density channel. Moreover, a zero sound mode 
appears in the spin channel, consistently with the results of Ref. [61]. Comparing 
the right and left panels of fig. 5.3, it is apparent that the collective mode strongly 
depletes the incoherent response, and largely contributes to the response in the 
|q| ^ limit. Figure 5.4 shows the response for different values of the density. Note 
that all curves are normalized so that the maximum of the corresponding Fermi gas 
structure functions be unity. 
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Figure 5.3. Left panel: spin-density structure function of neutron matter at p = po = 
0.16 fm~^, as function of A = cj/(|q|vi?. The curves have been obtained including only 
the Landau parameters with ^ = 0, 1 or taking into account the contribution associated 
with ^ = 2. For comparison, the thick dashed line shows the results of Ref. [61]. The 
curves are normalized as discussed in the text. Right panel: same as in the left panel, 
but for the density structure function. 




Figure 5.4. Left panel: spin-density structure function of neutron matter, as a function of 
A = cj/(|q|vi?, at different densities. The calculations have been carried out including 
the Landau parameters with ^ = 0, 1 and 2. For comparison, the thick dashed line 
shows the results of Ref. [61], corresponding to p = po = 0-16 fm~^. The curves are 
normalized as discussed in the text. Right panel: same as in the left panel, but for the 
density structure function. 
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5.4 Sum Rules in the |q| limit 

In the extreme low-momentum hmit, the dynamic response function is heavily 
affected by the static properties of matter. Linear response theory provides the 
qualitatively behavior of the the matrix elements of both the density and spin-density 
operator, as they must reproduce the static observables. Moreover, some conservation 
laws can be rearranged in such a way as to constrain the response function. The 
number of non-trivial sum rules grows with to the complexity of the hamiltonian, but 
they are usually limited by the degree of accuracy of the approximations employed 
in the calculations. Typically, these constraints are expressed by moments of the 
correlation functions, defined as 

— lim / Ljj'^Y(a,(jj)d(jj , 

and commutator relations. 

Consider a hamiltonian including central interactions only. The corresponding 
response function turns out to be diagonal in both the spin and density channels. 
The value of the moments with n odd can be simply expressed, e.g. in the spin 
channel, as 

...") (5.30) 

n times 

where the value of C, given by C = 27r/p in our case, depends on the pre- factor in 
the definition of the correlation function [62]. The first quantization form of cr is 

CT(c?) = ^crje''»-J . 

J 

As stated in Ref. [62], Landau theory does not exactly satisfy the spin-density sum 
rule for moments with n > 3. Moreover, in our case we the hamiltonian includes 
tensor interactions, although the corresponding Landau parameters are neglected. 

An alternative consistency test directly relates both the compressibility and the 
susceptibility to the corresponding dynamic form factors 

Our strategy is based on evaluating the left hand side of the above equations, 
separating the contributions of the collective mode from the one arising from single 
ph excitations, and then comparing it to the right hand side. 

In fig. 5.5, we show the behavior of ^~/^cr' "^^^ solid lines represents the right 
hand side of eq. (5.31) and (5.32), while the dots correspond to the left hand side. 
The results have been obtained for a neutron matter density p = 0.16 fm~^, for 
both the density and spin density channels. We can clearly see the two different 
contributions coming from single quasiparticle and collective mode excitation. As 
already pointed out, the latter gives a large contribution in the static limit. 
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(5.32) 
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Figure 5.5. Left panel: The sum rule in the density channel. Right panel: Same as in 
the left panel, but for the spin channel. In this case, the contributions of the collective 
mode (dotted line) and incoherent excitations (dashed line) are displayed separately. 
The observed trend reflects the one shown in fig. 4.3, showing that the collective mode 
plays a relevant role in the low momentum limit. 



Let us now consider the n — 1 sum rules. The integral of the Landau theory 
response in the left hand side can be evaluated explicitly both for the density and 
spin-density channels (see fig. (5.6)) . Within this approach we get, 

/2 

^^Landau(q'^)'^'^ = (5-33) 
m 

ml = |lim I uxi:nd.u{^,^)duJ = ^(1 + ^G^) • (5.34) 

While in the density channel the above expression is exactly the same as the one 
obtained from eq. (5.30), the corresponding expression for the spin channel turns out 
to be slightly different. The right hand side is still depending on both the symmetric 
and antisymmetric Landau parameters, and is no longer a function of the bare mass 
and the momentum transfer only (vrq^/m). This is intrinsically related to the degree 
of approximation inherent in Landau theory itself. In [44], the difference between 
eq. (5.30) and the values resulting from Monte Carlo simulations has been used 
to estimate the relative importance of multiparticle excitations in the kinematical 
regime where cj > |q|. 
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q (fm-^) 

Figure 5.6. The ^^/^ sum rule at p = 0.16 fm~^, as a function of momentum transfer. 
The sohd hne represents the right hand side of eqs. (5.33) and (5.34), while the dots 
correspond to the left hand side. The observed behavior do not change significantly with 
the density. 
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Chapter 6 

Neutrino mean free path 



The neutrino cross sections in nuclear matter plays a central role in determining 
neutrino transport in both supernovae and neutron stars. Many different scattering 
and absorption processes can occur and , at least in principle, they should all be 
taken into account. However, the importance of a specific process mainly depends 
on two distinct factors: the energy of the incoming neutrino and the composition of 
nuclear matter, which in turns depends on the density regime we are interested in. 

A list of the potentially relevant neutrino processes can be found in several 
papers, see, e.g., Refs. [64, 65, 66]. They include neutral current neutrino scattering 
off nucleons, alpha particles and nuclei, first discussed in the pioneering work of 
Freedman [67], super-allowed charged-current neutrino and antineutrino absorption 
on nucleons, neutrino-neutrino scattering and neutrino-antineutrino absorption. 
The inverse of several neutrino production processes, such as bremsstrahlung of 
neutrino-antineutrino pairs and direct and modified URCA processes, also contribute. 

An incomplete list of the processes that may be relevant at relatively high density 
includes 

z/g + ^ ^ e~ + p , z/g + P ^ e+ + n , z/^ + A ^ z/^ + A* , 

U£ + n ^ U£ -\- n , UQ + n-\-n^i~+p + n, U£ + A ^ i~ -\- A"" . 

It should be noted, however, that in the density region corresponding to the neutron 
star core, where the ground state of the matter is likely to be a neutron liquid, the 
occurrence of many processes may be inhibited by the absence of nuclei, or of a 
significant proton fraction. Typically, purely leptonic processes, such as neutrino- 
lepton scattering, are negligible in neutron stars. It has been shown that for neutrino 
with energy E^, ^10 MeV the charged-current cross section is two orders of magnitude 
larger than the z^e — e scattering cross section [65]. Nevertheless, this process may 
be important for the thermalization of emergent muon neutrinos [68]. 

Changing the initial conditions gives rise to several different processes. For 
example, at low temperature (T < 3 — 5 MeV) and relatively low density {p 2^ 
10^^ — 10^^ g/cm^) heavy nuclei are expected to be present, and dominate the neutrino 
cross section due to coherent scattering. In these conditions, the interparticle distance 
is large d 2^ 20 — 40 fm compared to the range of nuclear forces, the contribution 
of which is small, and correlations between particles are driven by the Coulomb 
interaction Fq — Z'^e^/d^ with Z 25. For low-energy neutrinos these coherent 
scattering processes are far more important than those involving nucleons [67]. 
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As the density increases, the nuclei get larger, and their separation distance 
descreases. Novel heterogeneous phases of matter, dubbed "pasta phases", have 
been predicted to occur: the nuclear surface begins to provide a significant con- 
tribution to the energy, and nuclei deform progressively from spherical to rod-like 
and slab-like shapes, in order t arrange themselves in a more energetically favorable 
configuration [69]. The occurrence of the pasta phase has been recently confirmed 
by a quantum molecular-dynamic inspired model [70]. A further increase of the 
density leads to the appearance of a homogeneous nucleon liquid. In this regime, 
low-energy neutrino interactions are strongly affected by nuclear matter dynamics, 
leading to the occurrence of many-body excitations, such as resonances and collective 
modes [71]. This work is focused on this density region. In the following Sections we 
will discuss the interaction of low-energy neutrino with neutron matter. Since the 
timescale of nuclear dynamic is much shorter than that associated to weak processes, 
we will assume that matter be in equilibrium and describe many-body effects within 
the framework of linear response theory using the effective interaction described in 
Chapter 2. 



6.1 Neutrino-neutron interactions 

The interactions of low-energy (typically few MeV) neutrinos can be described 
within the low-energy limit of the standard model of elect roweak interactions. The 
corresponding charged- and neutral-current interaction lagrangian densities read 

C'^- = ^l^j!^ for v, + n^i + X 

>C"^ = ^Wz for Vi^n^vt^X, 

where Gp — 1.436 x 10~^^ erg cm~^ is the Fermi coupling constant. While the 
leptonic current is well defined in terms of the degrees of freedom of the fundamental 
theory, the hadronic current involves composite particles and its form depends 
on the four momentum transfer associated with the process. However, when the 
latter is very small compared to the nucleon mass of ^ the hadronic current can 
be approximated with the simple form first proposed by Fermi to explain neutron 
/3-decay. Hence, we can write 

V = ^(1 - 75)7^^^ , iw = ^(^)(^v - 9Pah)l^n{x) 
Ifi = - 75)7^^^ , Jz = X{x){cy - CAj5hfin{x) , 

where X labels the final hadronic state and cy, gv? and gA denote the vector 
and axial coupling constant 



6.1.1 Neutrino- nucleon cross section in vacuum 

Lets us start by analyzing in detail the inclusive cross section of only process we 
need to consider in pure neutron matter 

iyi{k) + n{p) vi{k') + n{p) . 

^And also smaller then 1/L, L being the length scale specifying the size of the target particle, as 
pointed out in Ref. [65]. 
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We will first focus on the description of the above reaction in free space, and 
include medium modifications at a later stage. At leading order in perturbation 
theory, the invariant amplitude reads 

M = ^n,(k', 5)7^(1 -75K(k,s)(n(y)|J^|n(p)), (6.1) 

the four momenta of the participating particles being defined as 

p = iE,p),k^ (e,k),y = iE',p'),k' = (e',k') . 

Using standard techniques (see, e.g. Ref. [72]) one can write the differential cross 
section in the Lab frame in the form 



da 



m-p) 



where M denotes the neutron mass and the bar over refers to fact that the 
squared transition amplitude is averaged over the spins of the initial state particles, 
s and a and summed over the spins of the final state particles, and (J^ 
Substitution of eq. (6.1) into eq. (6.2) yields the expression 



that can be rearranged in the concise form 

£a{e,p) _ Gl 2Me' 



-L'^'^W^., (6.3) 



de'dQ (27r)2 16p-p)2]i/2^ 
where L^'^ and W^i, are the leptonic and hadronic tensor, respectively, defined as 
L^^ = Slk'^'k" + k^k'" - g^"'{k' ■ k) - ie^'"'^k'^k,3], (6.4) 

and 



6.1.2 Medium effects 

Non relativistic Fermi gas model 

We now want to improve upon the differential cross section (6.3), describing neutral 
current neutrino-nucleon scattering in vacuum, by including the effects of matter. In 
the region of low momentum transfer we can rely on the non relativistic approximation 
for nucleons, implying that the term arising from the fiux simplifies to [(k • p)'^]^^'^ 
Me 
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Within the Fermi gas model, neutron matter is seen as a collection of non 
interacting neutrons, the momentum of which is distributed according to the Fermi 
distribution. This model can be easily implemented in the formalism described in the 
previous Section, since the system can still be described in terms of single-particle 
states. 

The the hadronic current operator is applied to a state describing a neutron 
carrying momentum p, and the cross section is obtained collecting the contributions 
of all target particles, weighted with the Fermi distribution /(p) according to 

V ^ f 

or, more briefly 



de'd^ (27r)3 46 
the nuclear tensor being given by 

X (27r)5(£p+q - £p - w)5^{'k! + - k - p), (6.5) 

where p = N/V and the factor [1 — /(p')] describes the phase space above the Fermi 
level, available to the final state neutron. Rewriting the above expression in terms 
of the momentum and energy transfer, q = k — and w = e — we find 

= ^ ^ / (^-^^P^t^ " -^^P + ^)]Hp)\rzHp + q)){n{p + qMHp)) 
X (27r)(5(£p+q -Ep-U) . (6.6) 

In the non relativistic limit, 

and the only surviving components of the matter tensor are the density-density (00) 
and spin-spin (ii) ones 

Woo = Sid, CO) , T^^^" = -S(q,a;)% W^^'' = , 

while all off diagonal contributions vanish. Let us focus on the density channel first. 
Defining 

Pp(q) = (p|^'^^|p-q), 

the expression of ^(q,^^) can be written in the form 



= J/^e-*(p(q,t)p(-q,0))^G, (6.7) 
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where p(q, is time-dependent operator in the Heisenberg picture and 
denotes the expectation value in the Fermi gas ground state. 

In the spin-spin channel we find a similar result (note that as all diagonal terms 
are the same, we can choose z = 3) 

S\cL,u) = -J^ e-*{c73(q,t)c73(-q,0))re • 
Including interactions 

Once we have obtained the expression of the cross section in term of a ground state 
expectation value, including the interactions is formally easy. One has to make the 
replacement {. . .)fg with (. . .)int^ meaning that the expectation value is now in the 
ground state of the interacting system. We remind the reader that we are allowed 
to do so since neutrino physics and nuclear dynamics correspond to different time 
scale, and can therefore be decoupled. The length scale of the perturbation induced 
by the neutrino is much larger than the mean free path of a neutron in matter. As 
a result, the neutrino interacts with the system in equilibrium, and the zero sound 
mode can be excited. 

Contracting of the lepton and matter tensors yields 



(1 + COS 6>)S' + (3- COS e^)^^ 



where 9 is the angle between the directions of the incoming and outgoing neutrino. 
The resulting expression of the differential cross section is 

and the neutrino mean free path is given by 
^ = (6.8) 
GIp 



- /d^k'[(l + cos^)5(k-k^e-eO + (3-cos^)5A(k-k^e-eO] . 

4(27r)'^ J 

Note that, from a purely formal point of view, the above expression can be read- 
ily generalized to the case for hot neutron matter, substituting the ground state 
expectation value entering the definitions of S and Sa with an ensemble average. 

The nuclear correlation function S and ^Sa are related to response function 
through the fluctuation-dissipation theorem [49]. Consider an external probe P{-x^t) 
coupling to an observable of the system, say the density, through / d:s.p{-K,t)P{-K,t) 
then the following relation between the correlation S and the response function K 
holds: 2 ^ 

S{ci,uj) = ^Imi^(q,a;) . 



It follows that we can obtain the mean free path at T 7^ from 

+(3 - cos 6')Im[x''''(k -k',e- e')] }• (6.9) 
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Phace space considerations 

In this Section we review briefly outhne the calculation of the phase space integrations 
[61]. In an isotropic system the correlation functions only depend on the magnitude 
of the momentum transfer momentum q— |k — k^| and the energy transfer uj — e — e^ 
Choosing the 2;-axis along the direction of the momentum of the incoming neutrino, 
k, 9 of eq. (6.9) is the angle between k and k^, which can be written as a function 
of q and uj according to 



q^ = ^ k'^ - 2kk' cos 6 — ^ cos ( 



1 



+ 



2k'^ 2 /c2 



(6.10) 



In eq. (6.9), after the trivial integration over the azimuthal angle we are left with 
27T q^dqd cos 9 with 

cos 6qk 



1 



2k ~ 2k~q ^ 'q 



as required by energy and momentum conservation. Changing variables from (g, cos 9) 
to (g, uj) and taking into account the Jacobian of the transformation we find 

dqdcos9^ — \J{ ^^^^^^ )\dqdw 



dqw 



with 



dw 



(1- 



w 1 



(6.11) 




Figure 6.1. Left panel: Qualitative sketch of the phase space available at T = (green 
region). The light gray area represents the region determined by the degenercy condition 
only, eq. (6.16). Right panel: Same as the left panel, but for T > 0. It clearly appears 
that the available extends into the cj < region, corresponding to de-excitation of the 
system. 



In conclusion, we obtain an expression in terms of incoming neutrino momentum, 
fc, momentum transfer, and energy transfer uj 



1 



f dqduj / 
] (1 - e-P^) \ 



l{k) 2(27r)2 J {l-e-^^) 

X [(1 + cos e{q, UJ, k))K{q, oj) + {3 - cos 9{q, uj, k))x'"'{q, oj)] , 



(6.12) 
(6.13) 
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where cos 0{q, k) is given by eq. (6.10). The integration range is dictated by energy 
and momentum conservation. In the case of the incoherent excitation spectrum we 
have 

q = p^-p = k-k\ (6.15) 

Substitution of obtained from eq. (6.14) into eq. (6.15) allows one to determine 
the limits of the integration over the energy transfer cj 

-qv^i^<^<qv^^— , ^ = — • (6.16) 
2m* 2m* m* 

For a degenerate system, the maximum value of v is the Fermi velocity, vy-, and the 
term g^/2m* is usually neglected, with respect to qv^- A second constraint on cj, 
arising from neutrino kinematics, is 

\uj\<q<\uj- 2k\ . (6.17) 

Equations (6. 16)- (6. 17) determine the integration region for incoherent single-particle 
excitations. 

In the spin channel we also have the contribution of collective modes. Here the 
integration range is determined by eq. (6.17) only, since the dispersion relation of 
these excitations is fixed: w — -^vcq^ where zb corresponds to phonon emission or 
absorption. The resulting integration limits are 



2k 

0<q< — , uj = ±vcq 

l±vc 



Note that, so far, neither degeneracy nor temperature have been taken into 
account. Let us first include degeneracy at T = 0. Now all momenta below are 
filled, and there are no quasiparticles above the Fermi level. Although the system 
can be excited, with cj > 0, de-excitation is suppressed by Pauli principle. In this 
case we should add the condition (see fig. 6.1) : 

> . (6.18) 

Turning on the temperature, the energy transfer is bounded from below, as de- 
excitation is possible only in a range of energies ~ kT around the Fermi energy. As 
a consequence, 

-kT <uj . (6.19) 

The integration region is determined by the intersection of the above conditions. 

Figure 6.1 qualitatively illustrates the integration regions for both cold and hot 
neutron matter. A non zero temperature results in the extension of the available 
space to the cj < region, corresponding to de-excitation of the system. As the 
temperature increases, the lower limit can change appreciably, depending on whether 
or not the condition -kT > -qvp+q^ /2rrf is fulfilled. When -kT < -qvp+q^ /2rrf , 
degeneracy can also play an important role in limiting the phase space. 
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6.2 Numerical results 

This Section is devoted to the discussion of the numerical results of our analysis of 
the neutrino mean-free-path, A. 

Figure 6.2 shows the density dependence of the mean free path of a non degenerate 
neutrino with an energy e = 1 MeV in neutron matter at T = 0. The results have 
been obtained from Eq.(6.8), using the density and spin-density structure functions 
computed using the Landau parameters and Gi listed in Table 4.1 with ^ = 0, 1 
(solid line) and ^ = 0, 1 and 2 (thick dashed line). Comparison with the mean free 
path in a free neutron gas, displayed by the dashed line, shows that inclusion of 
interaction effects leads to a large enhancement of A over the whole density range. 
This behavior is to be ascribed to short range correlations arising from the repulsive 
core of the NN interaction, that prevent two interacting nucleons from being to 
close to one another. 




0.05 0.10 0.15 0.20 0.25 0.30 

p(fm-3) 



Figure 6.2. Density dependence of the mean free path of a non degenerate neutrino with 
an energy e = 1 MeV in neutron matter at T = 0. The density and spin-density structure 
functions have been computed using the Landau parameters and of Table 4.1 with 
^ = 0, 1 (open dots) and ^ = 0, 1 and 2 (solid line). For comparison the dot-dash line 
shows the mean free path in the free neutron gas. 

The left panel of Fig. 6.3 shows the energy dependence of the mean free path in 
neutron matter at density p — 0.16, corresponding to the Fermi temperatures Tp = 35 
MeV. The upper and lower curves have been obtained setting the temperature to 
T = and 2 MeV, respectively. The right panel also displays the density dependence, 
in the range 0.04 < p < 0.32 fm~^. It appears that the cross section exhibits a a ^ 
dependence on the energy of the incoming neutrino, and that at approximately e ^ T 
thermal effects begin to be less important. At e » T the curves corresponding to 
T = and T — 2 MeV are very close to each other, implying that the temperature 
dependence is neghgible. Moreover, due to the stiffness of the EOS, thermal effects 
appear to be less pronounced in the high-density region. 

The dependence of the mean free path of a non degenerate neutrino with energy 
e — 1 MeV upon both temperature and matter density is illustrated in Fig. 6.4. 
Note that, as the the density range 0.04 < p < 0.32 fm~^ corresponds to Fermi 
temperatures 14 < Tp < 55 MeV, the coUisionless condition T << Tp is always 
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;to(MeV) ^o(MeV) 

Figure 6.3. Left panel: Energy dependence of the mean free path of a non degenerate 
neutrino in neutron matter at different temperatures. Right panel : Same as in the left 
panel, but for different neutron matter densities. 



satisfied. We note again that in the high density region the EOS is very stiff and 
thermal effects become important at larger values of T. However, it must be pointed 
out that the pattern emerging from Fig. 6.4 results from the combined effects of 
several different factors, such as stiffness, phase space and degeneracy. 




lOOoL 



0.05 0.10 0.15 0.20 0.25 0.30 0.05 0.10 0.15 0.20 0.25 0.30 

p(fm-3) p(fm-3) 

Figure 6.4. Density dependence of the mean free path of a non degenerate neutrino with 
energy e = 1 MeV. The curves are labelled according to the values of temperature. The 
open dots in the left panel correspond to T = 0.01 MeV. 

Nevertheless the rate of depletion depends on density as shown in fig (6.4). In 
the high density region, the EOS is very stiff and particles are very packed one each 
other. As a consequence the effect of the temperature will start at higher value of T 
while we expect a deeper effect at low density. Several different factors like stiffness, 
phase space, degeneracy etc, inter-play at this level for T < and the trend is shown 
in the left panel of (6.4). 

As the temperature increase T > the available phase space saturates both for 
low and high density. The effect of the temperature depend now only on the factor 
in eq. (6.9), no sharp change are possible and A turns again to decrease with density. 
The right panel in fig (6.4) indicates the trend at high temperature. 



69 



Chapter 7 

Three-nucleon forces 



Over the past two decades, thanks to the availabihty of interaction models based 
on effective field theories, combined with the progress of many-body techniques, 
striking advances have been made towards the solution of nuclear Schroedinger 
equation. Several ab initio studies have been successfully carried out for light nuclei. 
Green's function Monte Carlo and No Core Shell Model calculations, performed using 
realistic potentials, provided exact solutions for A < 12 . For heavier nuclei, methods 
like the Self Consistent Green Function (SCGF)^, Coupled Cluster or sn-medium 
Similarity Renormalization Group (SRG) yield results that are still restricted to 
closed shell plus the addition/removal of one or two nucleons at most. G-matrix 
perturbation theory and the variational approaches based of CBF, such as the Fermi 
Hypher-Netted Chain (FHNC) scheme, have been extensively applied to infinite 
nuclear matter. 

The generalization of the available theoretical approaches to describe open shell 
nuclei, take into account the continuum and use realistic hamiltonians, including 
three-nucleon forces (3NF) is now regarded as a most prominent issue in nuclear 
theory. Recently, a technique to successfully attack the first two problems has been 
developed in Ref . [75] . This Chapter is devoted to the the discussion of the third 
one. We will provide a brief outline of the effects of 3NF illustrated in Fig. 1.7, and 
focus on how their inclusion affects the theoretical neutron drip line [10, 76]. 

Three nucleon interactions are a needed element of any refined calculations of 
both atomic nuclei and infinite nuclear matter. Phenomenological potential models, 
which are the most widely used in literature, allow to explain several important 
properties, from saturation of isospin symmetric matter to several features of nuclear 
spectroscopy. A very effective procedure to include the effects of 3NF is based on 
a density dependent modification of the two-body potential, inspired by the well 
known Fujita-Mijazawa two pion exchange process. 

In most of this Thesis work, we have followed the above scheme, as described in 
Section 2.3, taking into account both the attractive and the repulsive contributions 
to the 3NF in the definition of the effective interaction in a largely phenomenological 
fashion. In this Chapter we discuss a more microscopically rooted approach, in 
which the density-dependent correction is obtained in a fully-consistent background 



^ Which is also appHcable to study spectral functions and optical potential in a broad range of 
energies. 
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provided by by the Green function formalism. In the second quantization picture, 
the hamiltonian driving the dynamical interaction can be expanded in terms of 
creation/annihilation operators according to 

H = T + V + W = Ta^ala(3 + \ Va(3,jd at o^l asa^ + ^ Wa^j^seC ctl o^l ci\ 

t x \ I f - - + - - 4 

where T , V and W represent the kinetic energy and the bare two-body and three- 
body interactions, respectively. All the indexes run over the single particle states 
which are created/destroyed, depending on whether they are associated to a^a. The 
coefficients are completely defined imposing the full anti-symmetry of the matrices 
y, W. Note that the description in terms of creation/annihilation operators is the 
framework in which the Green function formalism has been originally developed. 
Here we use a state-of-the-art potential derived through a similarity renormalization 
inspired transformation (SRG) [77]. The SRG evolution is an unitary rotation that 
allows to transform a given hamiltonian into a softer, low-momentum, interaction, 
suitable for use in a perturbative approach. Two different initial hamiltonians based 
on xPT have been used: a two-body potential calculated up to N^LO and a N^LO 
two-body plus a N^LO three-body potential. The flow parameter characterizing the 
evolution has been fixed to o: = 0.8 fm~^ and we have retained all the evolved terms 
up to three-body level, getting an induced (labelled with ind) 2 + 3-body force in 
the former case and a full (labelled with full) 2 + 3-body force in the latter [78]. 

7.1 Green function formalism 

Green's function techniques are well established as a very powerful tool to describe 
quantum systems in a variety of fields of Physics. This formalism proved very 
effective both in vacuum and in matter, as well as in and out equilibrium. In 
this section, the Self Consistent Green Function (SCGF) approach is applied to 
many-body theory. In this context, it is capable to provide a clear description of 
the dynamics, including a wide range of observables, such as binding energies and 
spectral functions, which are analyzed using a diagrammatic technique. In the case 
of non-perturbative dynamics, a self-consistent approach is still possible within the 
equation of motion method. The fundamental quantity of this formalism is the 
single particle Green function, that describes the propagation of a particle, or a hole, 
within the system 

9a^it - t') = -i«|rK(t)a^(t')]K) . (7.1) 

In the above equation, is the ground-state wave function, T is the time-ordering 
operator and aa/al^ denote the particle annihilation/creation operator in the state 
a. The quantity defined by eq. (7.1) describes to the creation of a single particle 
state /3 at time together with the destruction of a single particle state a at time t 
for t — < 0, and the destruction of a state a at time t together with the creation 



7.1 Green function formalism 



71 



of a single particle state /? at t' for t — t' < {). Thanks to energy conservation, 
this function depends on {t — t') only. It is remarkable that once this function 
is known, all one-body can be described^. This can be easily seen in the second 
quantization formalism, in which a one-body operator can be expanded in terms 
of creation/annihilation, a^o) ^ operators according to O^^ — ^C6(3 ^1x^(3- ^ 
results, the expectation value is given by 



A clear picture of the information contained in is provided by the the Lehmann 
representation, obtained transforming to Fourier space and using the completeness 
relation fulfilled by the eigenstates of the N zb 1-body system: 



u - {Ek'^' -E^) + if] ' u - « - E^^-') - if] 



The poles of the above equation provide for the excitation spectra ±{E^^'^ — 
Eq^) asociated with addition/removal of a particle to/from the ground state. The 
corresponding residues, reflecting the transition amplitudes, go under the name of 
spectroscopic amplitudes. The hole part of the propagator gives information on the 
process of particle emission, the poles being the exact energy absorbed in the process. 
Note that this analytic structure is completely general and that these energies and 
amplitudes are obtained solving a Schroedinger-like equation, that in diagrammatic 
language is usually referred to as Dyson equation 



K^) = ^LK^) + 9a|3{^)^|3'{^)9|3^a^{(^) • (7-3) 



In the above equation, is the unperturbed, free propagator, corresponding 
to nucleons moving without experiencing "dynamical" interactions. However, 
statistical correlations, indiced by Pauli exclusion principle, are always present. 
The matrix Ti'^{(jj) embodies all the dynamical information stored in the Green 
function, and the star refer to the fact that only irreducible diagrams, i.e. dia- 
grams that cannot be obtained combining lower-order diagrams already present in 
S'*'(cj), must be included. Equation (7.3) can be written in diagrammatic form as 



^This is also true for some peculiar multi-particle operators. 
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Diagrammatic representation of Dyson equation (7.3) (left), and its perturbative 
expansion (right). 

It is wort pointing out that the full knowledge of would be equivalent to 

solve the Schroedinger equation without any approximations. In actual calculations, 
we need to select a restricted number of irreducible diagrams, which amounts to 
limiting many-body correction to a maximum number of particles. As already 
stated, dynamical effects come from the self-energy S*, which, in principle, includes 
all the diagrams that can be seen as a one-body correction arising from the bare 
hamiltonian. This is what has been usually done with hamiltonians which include 
two-body interactions only. Switching on three-body interactions, the convergence 
of the Dyson equation becomes very slow, since the thee-body term is actually 
transformed into a two-body effective interaction. The corresponding diagrams can 
appear at any order of the expansion, and heavily affect the results. However, thanks 
to the implicit form of the Dyson equation one can find an algorithm to sum up 
these diagrams at all orders. The basic idea is rearranging the diagrams with the 
help of an effective hamiltonian i7, defined as 



H = Ta(3 + i Va(3,jS da ^(3 CiS + ^ ^afSj.SeC Cia ^(3 ^7 ^'e ^5 



f-- + --4 



The new hamiltonian has a clear physical interpretation: T, V are effective operators 
which include the hare T and V terms as well as the screening arising from the 
three-body force. Note that H has the same dynamical content as i7, with the 
constraint that now W must act only as a pure three-body interaction, in order to 
avoid double counting of the screening effect. In Ref. [76] it has been shown that the 
dominant effect of the three-body force can is in fact included using T and while 
the pure three-body term plays a minor role. In the following, we will disregard the 
last term, the importance of which has been critically reviewed, up to third order, 
in Ref. [81]. 

In diagrammatic language the two operators T and V take a transparent form, 
in which gPP/^^ is the two-body propagator, that can be defined through a straight- 
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forward extension of eq. (7.1), and the three-body propagator also appears. 





j|rv\/s/\/N/v/\^ ^ — |— — — — 




Figure 7.1. Diagrammatic representation of the one and two-body contributions to the 
effective hamiltonian H. 



gika'P'A^ - t') = -i{*^|rK(t)a^(t)a«(t)al,(t')al;,(t')4'(*')]l*^) (7-4) 

Proceeding as in the derivation of eq. (7.2), we can define the average of two- and 
three-body operators 

\^ / — ^a(5,a'^' Pa'(3',a(3 Pa^a'(3' " \^0 /^^l / 

{^U ) — ^a(3-i,a'(3'-f'Pa'l3'-i',al3-i Pa(3-fa^(3^Y \^0 \(^a'^l3'%'^l^^^(^\^ / ^ 



with 



The diagrammatic representations of the corresponding expectation values are 



(o 



1B\ 






Figure 7.2. Diagrammatic representation of the expectation values of many-body operators 
evaluated through the full Green functions. 



7.2 Sum rule for the binding energy 

In principle, the calculation of the binding energy of a nucleus binding energy within 
the Green function formalism is not an easy task . The evaluation of {H) directly 
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from kinetic, T, and potential, V and W terms, involves operators acting on states 
describing an up to three particles. As shown in fig. (7.2), one would have to 
calculate independently three different expectation values, involving three different 
propagators from three different Schroedigner-like equations. This would of course 
imply a large amount of computer time, as well as a huge theoretical effort. 

In this Section, we will show how the number of Green functions involved in the 
calculation can be reduced, exploiting a theoretical result first obtained by Koltun. 
The underlying idea stems from the observation that the operator we want to average 
is the same which drives the time evolution of the system. In the Heisenberg picture, 
we can write the time evolution equation for the annihilation operator in the form 



(7.5) 



Defining the time-derivative of (7.1) as 



dga/3{t) 



dt 



= i(*ol«*(0)^*^ 



l*o) 



and using the relations 



ZTTl Jc^ 

^at^[a,,r] =r, 

Y,al[a^,V]^2V, 

OL 



(7.6) 



(7.7) 



we obtain 



(r + 2^ + 31^) =^(*o|aiK,ii-]|*o) 



:E(*o|al(o) 



da, 



exit) 



dt 



t=o+ 



^ 2tii Jc^ 



(7i 



Equations (7.5) and (7.7) lead to eq. (7.8), which is the extended Kultun sum rule, 
modified from the original result to include three nucleon forces. The C \ integration 
contour must be closed the in the upper half-plane, in order to extract the residue 
of quasi-hole pole. 

Different extrapolations of the ground state energy {H) = {T + V + W) can be 
inferred from the above equations. The one we use is 

(^) = E^/^ d^[Tap + ^S^p]9ap{uj)-\{W) , (7.9) 



a/3 
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where the average of W must be calculated separately. A similar formula, involving 
thetwo-body operator V , can also be obtained 



In the following, we mainly refer to eq. (7.9) since {W) / {V) ^ 10%. Consequently 
the contribution of (W) would have a minor impact in the final energy. 



7.3 Self-energy and iterative method 

The degree of accuracy of the SCGF formalism mainly depends on the number of 
irreducible diagrams one is able to include (and to sum up) in the self-energy. The 
approach we are going to described has been developed for infinite nuclear matter, 
and has been described in several papers for the caseof a hamiltonian involving a 
two-body interactions only [79, 80]. Here we extend the method taking into account 
three-body effects using the effective hamiltonian H. 

The diagrams involving the purely three-body interaction W have been discarded. 
As already stated, they should play a minor role [81]. We use a third order ap- 
proximation, referred to as Fadeev-Tamm Dancoff Approximation (Faddev-TDA) 
approach, in which the contribution of two particle-one hole {2plh) and two hole-one 
particle intermediate states are taken into account self-consistently, and the inter- 
action between pp/hh or ph excitations are modeled within the TDA scheme. The 
equations describing the polarization and particle-particle (hole-hole) propagator 
in TDA approximation are depicted by the yellow bubbles is fig. (7.3), where the 
effective V is employed. These polarization and pp/hh propagators need to be added 
consistently to determine the pp-h^hh-p propagator R{uj)^ yielding the third order 
contribution. 

Its contribution is sketched in fig. (7.3), for the case of the pp-h channel. 
It appears that it includes the effect of ph and pp/hh motion, allowing for the 
interferences between them and giving at the same time the correct combinatorial 
factor to each diagrams, without the need of subtracting spurious terms. Moreover, 
in the same figure we show how the new terms arising in T and V can be summed 
up within the same scheme. For a exhausting review, see Ref. ([82]). 

In the SCGF approach the self energy matrix S"*" is expanded in term of the 
dressed propagator up to the required perturbative order. This means that the actual 
degrees of freedom are the excitation of the fully correlated system, and the effects 
of fragmentation are already included in the iteration scheme for self-consistency. In 
applications, we need a first approximation to the propagator, e.g. the propagator 
obtained within the Hartree-Fock approximation, to start the calculation. Solving 
for the first time the Dyson equation with this propagator we obtain a new g{(jj)^ 
which is then used to calculate again the self-energy. This procedure is iterated until 
convergence is reached. 




(7.10) 
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7. Three-nucleon forces 





|»\/\/v/\/\/v/\^ ^= ^ — |— ^ ~ ~ ~ ^ 




Figure 7.3. Diagrammatic representation of the self-energy in the Faddeev-TDA approx- 
imation. The 2hlp contribution is not represented. The yehow boxes correspond to 
the Green function in TDA approximation. In this case the double line represents the 
self-consistent Hartee-Fock propagator. 



7.4 Approximations and numerical results 

Including the screening term in T and V amounts to using one- and two-body density- 
dependent interactions. This procedure sharply increases the required computational 
effort, since these interactions must be rebuilt at each interaction. We adopt several 
approximations aimed at making the code more effective and, at the same time, at 
minimizing the loss of information. In T and (VF), we approximate the two-body 
and three body propagators with two and three single particle propagators, which 
are by far the most accurately computed quantities within our approach. This 
approximation implies that the interaction between the single particles described by 
the Green functions are negnected. Moreover, we stop the evaluation of T and V at 
the first iteraction. 




Figure 7.4. Diagrammatic representation of the approximation scheme employed to carry 
out for the three-body average. The double line represents the full propagator in 
Faddeev-TDA approximation. 

The results obtained using the above procedure are collected in Fig. 7.5, showing 
the binding energy of the main oxygen isotopes. Consider the behavior of the 
energy of ^^O, represented by thick dashed line. It appears that when we use 
the full three-body interaction, ^^O is unbound. On the other hand, when the 
induced three-nucleon force is employed, it seems to be bound, with a ground state 



7.4 Approximations and numerical results 
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energy slightly below that of O. This is to be ascribed to the fact that the latter 
interaction actually plays the role of a pure two-body interaction, in which the 
repulsive component of the three-nucleon force is disregarded. The overall picture 
clearly indicate that a proper treatment of the three-body force is essential to 
reproduce the oxygen neutron drip line. 




456789 10 11 456789 10 11 



eMax eMax 

Figure 7.5. Behavior of the binding energies of the main oxygen isotopes, as function of 
the energy of the mayor shell employed, eMax. The results shown in the left and right 
panels have been obtained using the full three-nucleon force, interaction and the induced 
interaction, respectively. The thin dashed lines indicate the experimental values. The 
thick dashed line corresponds to ^^O, which appears to be unbound in the left panel and 
bound in the right panel. The results have been obtained using a harmonic oscillator 
potential with Tluj = 20 MeV. 
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Summary &: Outlook 



Over the past years, several refined simulations of supernovae have appeared in the 
literature [83, 84]. The main problem of these simulations is the description of the 
final stage of the evolution, that is the gravitational decoupling between the core 
and the mantle leading to the birth of a proto-neutron star. Many ideas have been 
proposed and already tested to explain the failure to achieve this goal. All clues seem 
to point to the oversimplified description of neutrino driven mechanisms in matter^ 
[85] . Indeed, neutrino opacity of nuclear matter is far from being understood at 
fully quantitative level. 

The impact of nuclear correlations on free the Fermi gas picture appears to be 
impressive, resulting in the appearance of a complicated pattern of excitations, and 
conventional perturbative approaches are designed to fail when dealing with nuclear 
matter structure and dynamics. 

The aim of this Thesis has been the development of a realistic description of 
neutrino interactions with nuclear matter. We have restricted our analysis to the 
case of pure neutron matter, and evaluated the neutrino mean free path as function 
of both density and temperatures near T ^ 0. In addition, we have extended our 
analysis to cover a broad density range, up to twice nuclear saturation density po- 
In fact, there are indications that the depletion of the hard repulsive core of the 
NN interaction allows the liquid phase to persist well beyond saturation density po, 
moving the threshold of the transitions to more "exotic" phases to densities as high 
as three- four times pQ. 

The starting point of our work has been the two-nucleon interaction in vacuum. 
We employed a phenomenological potential based on the Argonne vis model, yielding 
the best available fits of NN scattering data. Medium effects have been included 
through the CBF formalism, exploiting the expression of the energy at two-body 
level in the cluster expansion. This procedure provides an effective interaction taking 
into account medium effects, such as screening of the repulsive core arising from NN 
correlations. The main advantages of this approach is that the effective potential 
can be used to carry out perturbative calculations in the Fermi gas basis, and in the 
p ^ limit reduces to bare Argonne potential by construction. 

The Standard Model of Particle Physics, in the low energy limit, provides an 
extremely accurate description of semi-leptonic weak processes through a vector 
(CP conserving) and an axial (CP violating) currents, responsible for the so called 
Fermi and Gamow- Teller transitions, respectively. The description of these two 
transitions can be formally extended to the case on neutrino interactions in matter. 



^Of course, there could also be "new physics" that we are completely unaware of. 
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7. Three-nucleon forces 



although the relevant mechanisms strongly depend on the energy scale. In the 
region of low momentum transfer, many-body effects play a crucial role, but the 
self-consistent evaluation of the dynamic form factors within CBF involves non 
trivial computational problems. The authors of Ref. [48] succeeded in carrying out 
a consistent CBF calculation of the response restricted to the density-density sector 
only. 

We have followed a different approach, based on Landau theory of normal Fermi 
liquids. In this case the CBF effective interaction has been employed to obtain the 
numerical values of the Landau parameters. 

Within Landau theory the ring diagrams, which are needed to describe the 
response at low momentum transfer treating both the coherent and incoherent 
contributions on equal footing, can be easily summed up. While the static properties 
of neutron matter obtained in the Landau scheme have the same degree of accuracy 
as the previously available CBF results, within our approach we can clearly see 
the emergence of the zero-sound contribution. Moreover, our results are in fairly 
good agreement with the results of highly refined many-body calculation based on 
similar dynamic models. Our response functions also appear to be consistent with 
the existing literature. 

The evaluation of neutrino mean free path was the ultimate goal of the Thesis. 
We find that the dependence on neutrino energy, e, is roughly e^, consistent with 
available estimates and different from the dependence of the free Fermi gas model. 
Our analysis also included the dependence on density and temperature, limited to 
the region T < 10 MeV, where the coUisionless approximation is expected to be 
applicable. 
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Appendix A 

Properties of the operators 

A.l algebra 

The six operators already introduced Chapter 1 [/, (cri -(72), S'i2(f )] ® [/, (ti •T2)] have 
the very important properties of forming an algebra with respect to multiplication, 
as they satisfy the relations 

k 

The values of K^j can be easily found exploiting the SU(2) algebra of Pauh matrices 
[39] . The fundamental representation is commonly defined choosing (73 in diagonal 
form, 

and the commutation and anti-commutation relations are 

{ai^aj} = 26ij. 

A. 2 Two-particle system 

In order to derive, and then solve, the Euler-Lagrange equations for the correlation 
functions, it is convenient to change basis, and use operators that partially decouple 
in the total spin-isospin space. Let us consider two particles labelled by their spins 
(81,82) = (cri,cr2)/2, respectively. The total spin states \S,Ms) are related to the 
single-particle spin states by 

|oo) = -^(in)-ut)) 
|io) = ^(in) + iit)) 

|i - 1) = I U) 

111) = I tt). (A.l) 
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A. Properties of the operators O 



In operatorial language we can write 

S2 = S? + S2 + 2(Si-S2), 

implying {S = 0, l|cri • cr2\S = 0, 1) = —3, 1, respectively. 

The projection operators on states of definite total spin are 



Ps^o - Pi = i^l^:^ (A.2) 



4 

and satisfy the standard properties of a projector 

where I is the identity. 

We can easily construct the spin-exchange operator P^- = P3 — Pi such that, 

P,\SMs) ^ i-f^'\SMs) . 

The same relation can be obviously used in isospin space, so that we can define a 
total spin-isospin exchange operator P^r^ = PaPr such that 

P^rlSMs^TMr) = i-f^^lSMs^TMr). 
A. 3 Tensor term 

Let us now turn to the tensor operator 512, defined as 

5^12 = 3(0-1 • r)(o-2 • f) - (o-i • 0-2), 

where r is a unit vector along the direction of the relative coordinate of particles 1 
and 2, while r = |r|. Using the properties of Pauli matrices, it can be easily show 
that 

SuicTi • (T2) (cTi • Cr2)S'i2 Si2- 

This relation implies that ^12 acts on 5 = 1 state only and annihilates 5 = state-*^. 
Moreover, 

8^2 = 6 - 25i2 + 2(o-i • 0-2) = 8 - 2Si2 , 

and 

3 r 

VS'12 - 0-1(0-2 • f) + 0-2(0-1 • f) - 2f (0-1 • f)(o-2 • f) 

r L 

^^Si2 = -^Si2 . (A.3) 

The above relations imply that for a generic function u{r) depending on the radial 
coordinate r 

(Vu) ■ {VS12) = (^) {(Ti ■ 0-2) • (VSu) = . 



^Recall that the S = state is an eigenstate of (cri • (T2) with eigenvalue —3 



A.4 Total spin-isospin representation 
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Further useful properties of the tensor operator are 

6 



36 



i(S 



(V5i2)^ = £^(8-5i2) 
[5'l2, (VS'12)] 
[Sl2,{VSi2)]V 

(V5i2)[5i2,V] =-(V5i2)' 



= 



(A.4) 



A.4 Total spin-isospin representation 

In this section we give the exphcit expressions of the matrices needed to transform 
from the basis of the to the total spin-isospin representation. Consider a generic 
operator x written in the form 



x= ^x^(r)0" 



n=l 



= Xc + Xr{Ti • T2) + Xcj{(Ti • (T2) + X^ricri ' Cr2)(Ti • T2) + XtSu + XtrSi2{Ti • T2) . 

In the representation appropriate to describe channels of fixed spin and isospin it 
reads 

x = Yl [^TS + SsiXtrSu P25+in2T+i, (A.5) 



TS 



with 

/ 1 -3 -3 9 \ 
11-3-3 

1-3 1-3 

\ 1 1 1 1 

implying 



\ 


( 


Xq 






( ^00 \ 


















X(j 










[ 


Xq-t 


) 







1 -3 
1 1 



Xt 

XtT 



[ XTS =Xe + (4r-3)x^ + (45-3K + (45-3)(4r3)x, 
I ^xt + (4r - ?>)xtT 

Finally, the inverse transformation is 





/ 1 


3 


3 


9 \ 




( ^ 




/ 


Xq 




1 


-1 


1 


-3 


3 














16 


-1 


-3 


1 


3 










X(j 






I 1 


-1 


-1 


1 ) 




^ xii ) 




v 


X(JT 


) 



1 

4 



1 3 
-1 1 



xt^ 
xti 



Xt^ 

xti 



Xt 

XtT 
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Appendix B 

Energy functional 



The energy per particle at two-body cluster level can be written (see Eqs.(2.6) and 
(2.8)) 

^\fl2,[tl+t2, /12] 



i<j 



+ fuvufu \ij -ji) , (B.l) 



with 

ti = -^^V? , + t2 = -1 V2 - , (B.2) 

2m m 4m 

where V acts on the relative coordinate r, while Vi^ acts on the center of mass 
coordinate R, defined as 

r = ri-r2 , R=^(ri + r2), (B.3) 

respectively. 

Using the static part of the interaction, both the correlation function /12 and 
the two-nucleon potential vu are written as 

/i2 = E/^(ri2)0?2 , t^i2 = Et'^(ri2)0?2, (B.4) 

p=l p=l 

with the six operator 05^2 the properties of which properties are discussed in Appendix 
A. 

The Fermi gas (FG) two-nucleon state is given by 

\ij) = ^e'^^^-''^^^^-''^^ \SMs,TMt) 

with 

|k^|, \kj\ < pf 



k=^(k,-k,) , K = k, + k, . (B.6) 



We will discuss the potential and kinetic energy term separately. 
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B. Energy functional 



B.l Potential energy 

Consider the operator 

m2 /l2'^12/l2 , 

and the decomposition of /12 in the TS-representation (see Eq.(A.5)) 



/12 = E 

ST 



fsT + SsiftrSu 



^2S'+in2r+i 



(B.7) 



In the above equation, 1127^^1 are spin and isospin projection operators, 

whose properties are given in Appendix A. Performing the decomposition for W12 
and V12 we find. 



^12 = ^YsofTo^To + ^si{vTi fri ^ ^ftT ^ ^ {friftT - ftr) S12 

in 



+ 



TS 

+ ^tT [JTlJtT - JtTj yJTl - "^JTlJtT "^^Jtl J ^12J / J ^2^+1^^2T+l ^ 

and then we can identify 

^Ti = ^Ti (/ti + + Wv^T [friftT - ftr) (B.9) 
= 2vj^i [fTiftT ~ ftr) + ^tri^fTi ~ ^friftT + ^'^fti) • 

1 



/to 



Now, after replacing 



i<j ij 

the potential energy contribution to {AE)2 reads 



(B.IO) 



1 1 

2^2 



E E EE 

SMs TMt k^k^- S-^T^ I 



55/i'?i;tT^(r)(S'i2P25^+in2r^+i) 



^5/T/(r)(P25^+in2T^+i) + 
Jd\id\2e'^^^"'-^^"'^ (B.ll) 



7/;^/T/(r)(P2^^+in2r^+iPar) + 5^a'^^tTK0(5'i2P2^^+in2r^+iPar) 



where Pcrr is the spin-isospin exchange operator defined in Appendix A and the 
expectation values (O) are taken over two-nucleon states of definite total spin and 
isospin \S Ms, T Mt). Using 



J dVidV2 = J d^rd^R = V J d^r , 
the definition of the Slater function. 



J2 e"^"^ 



V 



|k|<PF 



(27r)3 y|k|<p 



N 



(B.12) 



(B.13) 



B.2 Kinetic energy 
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and the results of Appendix A, we finally obtain 



21/2 



__y ^(25+ 1) (2r+ 1) / d\wsT{r) [l - {-lf+^£\p^r) 

^ ST ^ 



for example, in the case of symmetric nuclear matter (z/ = 4), 

+ [3woi{r) + 3^10 (r)]a+(pFr)] I 

where p = N/V is the density and 

a±{x) = l±f{x) . 

B.2 Kinetic energy 

Let us now discuss the kinetic contribution to the energy, given by 



(B.14) 



(B.15) 



(B.16) 



/l2,[tl+t2, fu] 



1 

2m 



/l2, V^, /l2 



(B.17) 



We consider spin-zero and spin-one channels separately. 



Spin-zero channels In these channels, the relevant part of the correlation function 
is given by 

/i2 = E/™W^in2T+i . (B.18) 

T 

Making use of the results of Appendix A, as well as the relation 



/to, vVto =0 , /to, (V/to)V =-(V/ 



we find 



/i2, v^ /i2 



ITO 



(B.19) 



TT' 

= ^ [/to, V^, /to Pi^2T+1^2T'+1 
TT' 

= E[/™' (VVto)+2(V/to)V 

T 

= 2^ [/to, (V/to)V 

T 

= -2^(v/To)'Pin2^+i . 



Pi n2T+i 



Pi n2T+i 



Finally, 



1 

2m 



/l2, /i2 



m 



E (WTo)'Pin 

2T+1 



(B.20) 



(B.21) 



88 



B. Energy functional 



Spin-one channels In these channels, the correlation function is given by 



fTi{r) + ftT{r)Si2 



-P3 n2T_|_i 



(B.22) 



Relying once more on Appendix A, we calculate 

E [V"' ifT'l + ftT'Sl2)P3U2T'+l] = E {[V'' fT'l] + [V, /ir"5l2]}P3n2T 



'+1 



= E {(V'/t'i) + 2(V/tT')V + (V^ftT'Su) + 2(V/iT'5i2)v}P3n2r.+i 
T' 

= E {(V'/t'i) + 2(V/iT')V + {V^ftT')Su + {V^Si2)ftT' 
T' 

+ 2(V/tr')(V5i2) + 25i2(V/tr')V + 2/tT'(V5i2) v}P3n2r'+i • (B.23) 

Hence, the commutator in Eq.(B.17) can be rewritten as 



fi2, V2, /: 



12 



= E[/t1 + /*t5i2,{---} 

= EK^ + 4^0^3n: 



-P3n2T+l 



'■2T+1 5 



with 



and 



/ti,{...} 



p(2) 



ftrSu, {■ ■ ■} 



(B.24) 
(B.25) 



{ . . . } = {(vVt'i) + 2(V/jr')V + (vVtT')5i2 + (v25i2)/tT' 

+ 2(V/iT')(V5i2) + 25i2(V/tr')V + 2/iT'(V5i2)v} . (B.26) 



We find 
and 



5^(1) 

T 



-2(V/ri) -2(V/ri)(V/iT)5i2 , 



(B.27) 



i^(2) _ 



ftrSu, 2(V/ri)vJ + [ftTSu, 25i2(V/ti)V 

+ [ftTSl2, 2/ti(V5i2)V = 

= -2(V/ri) (V/ir)5i2 - 2{VftTfsj2 + 
+2f^Si2, (V5i2)V 

= -2 (V/Ti)(V/tr)<Si2-2(V/tT)'(8-25i2) + 



+ 



2/ir 



^(L.S) + -^(8-5i2) 



(B.28) 



B.3 Explicit form of {AE)^ 
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Collecting all pieces togheter, we find for the spin-one channels 



1 

2m 



/i2, V^/: 



12 



m 



E (V/ti)' + {VfTl){VftTSl2) + 



+ (V/tT) (8-25i2) +/|r 
1 



§(L.S) + ^ (8-5,2) 



>P3 1127^+1 



E (V/t5) + ^51 [2(V/t5) (V/tT)^12 + 
+ (V/,t)'^12 + ftT^ (L . S) + ^ (8 - ^12)] |p25+in2T+l 



= "^[trsir) + Ssi 



TS 



ttT{r)Si2 + hT{r){L'S) 



}P2S+1^2T+l ^ (B.29) 



with 



Ut 

tbT 



m 

^ '(V/ti)' + 8(VAt)^ + ^//t 



m L 
1 

m 
1 36 2 



2 , 48 ,2 
2(VM)(V/tT) - 2{VftTY - ^ f'rr 



.2 6 



B.3 Explicit form of {AE)^ 

We can rewrite 



{AE)^ = J2{ij\Wi2\ij-ji) 

i<j 



with 



Wl2 = 



1 



m 



/i2, v^/] 



12 



+ /l2'yi2/l2 



where 



T5 

Wto 



WtT{r)Si2 + WbT{r){L-S) 



^P2S+1^2T+l 



m 
1 

m 



(V/Tl)' + 8(VAT)' + ^/i2T 



+ 



+^ri (/ri + 8/tT) + 16?^tT(/ri/tr - ftr) 



m 



6 



W^tT = - 2(V/Tl)(V/tT) - 2(V/tT) - ^ //t 



+ 



W^6T = 



1 36 



.2 -^tr • 



(B.30) 
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B. Energy functional 



Making use of the expression for the expectation values given in Appendix A, 
we finally obtain for symmetric matter (compare to Eq.(B.15)) 



N 



P_ 
32 



Woo{r) + 9Wnir) a-ippr) + 



+ 3Woi{r) + 3Ww{r) a+ippr) } . 



(B.31) 
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Appendix C 

Euler-Lagrange equations 



C.l Spin singlet channels: uncoupled equations 

In the spin-zero channels, the energy per particle of SNM, evaluated at two-body 
cluster level, reads (compare to Eqs.(B.15) and (B.21)) 



(A^)^ 



32 

^ (2r + l) 47r J r^dr\{f!rof + VTof^o aroiPFr) 



POO 

= const / dr F 
Jo 

where aTs{x) = 1 — {—)^^^f'{x) and 
^ /to 5 /to 



/to 5 /to 



(/to) + "^Tofro 



^TO 5 



(C.l) 



(C.2) 



with 



0TO Ty/a^ 



(C.3) 



The corresponding Euler-Lagrange (EL) equations for the unknown functions /to 
are given by 

d dF 



From 



we obtain 



dr dffro dfr 



TO 



OF 

OF 

Wo 
d dF 



2 VTo fro <Pto ' 
2 /to '/'to 5 



2 /to '/'to + 4 /;^o '/'to '/'to 



TO 



J TO YTO 



TO YTO 



vto / 



TO </^ro 



. 



(C.4) 



(C.5) 
(C.6) 



92 



C. Euler-Lagrange equations 



Introducing 



9to = /to 0TO , 



we can put Eq.(C.6) in the form 



If 



(C.7) 



(Ci 



Now we introduce a Lagrange multiplier, in order to fulfill the requirement 
The resulting equation is Eq.(4) of [73] 



If 

9to 



6to 



+ {vto + A) ^7^0 = 



(C.IO) 



to be integrated with the boundary conditions 

5ro|.=o = , (C.ll) 

9T0\r=d = (t^T0\r=d ■ (^-12) 

C.2 Spin triplet channels: coupled equations 

In the spin-one channels, the contribution to the energy is given by (see Eqs.(B.15) 
and (B.31)) 



P_ 
32 



(VM)' + 8(V/,^)' + S/'t 



(2r + i) I dV| 

const / dr F f^^, ftj,; /^i, flj, , 
JO L J 



(C.13) 



where 



/ti5 /tT5 /ti5 ftT 



' a2 ,2 



+ 8(/, 



' ^2 1,2 



48 



f2 1,2 



Tl) V^Tl o Utrj V^Tl + "2 /tTV^Tl 



+ 



+ 



(/tiAt-At)] ■ (C.14) 



In this case we have two coupled EL equations 

= 

= . 

Carrying out the derivativees as in the spin-zero channels and defining 

9T1 = /ti0ti , 9tT = V8ftT(pTi , (C.15) 



d OF 


dF 


dr df^^ 




d dF 


dF 


dr df'^rj. 


^ftT 



C.2 Spin triplet channels: coupled equations 
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we find 



9'It 



9ti - (fe + ^Ti) 9T1 - VSvtrgtT = 

^ + {vTi - ^^vtr) + ^ 9tT - V8vtT9Ti = . 
Finally, inclusion of the Lagrange multipliers needed to guarantee 

9Tl\r=di ^ ^Tl\r=di ' 

/I _ J,/ I 

ytT\r=d2 ~ ^Tl\r=d2 ' 

with, in general, di 7^ ^2, leads to (compare to Eq.(5) of [73]) 
^ + {vTi + Ai) gri - (^VSvtT + A2) gtr = 



^^1 



9^T 



(t>Tl 



+ {VTl - '2'VtT + Al) + ^ QtT- [VSVtT + A2) STi = , 

with the boundary conditions 

9Tl\r=0 = ' 

9Tl\r=di = (t)Tl\r=di ' 



and 



9tT\r=0 = ' 
9tT\r=d2 = • 



(C.16) 



(C.17) 
(C.18) 



(C.19) 



(C.20) 
(C.21) 



(C.22) 
(C.23) 
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Appendix D 

Tensor operator in momentum 
space 



In [47] the Fourier transform (FT) of all the eighteen operators of Argonne potential 
are performed exploiting the equivalence between q and — iVr- Here we use a 
different approach, in which the FT is obtained carrying out the integrations through 
the decomposition over the orbital wave functions. 

We will start from the Suir) operator as it appears in the expression of the 
energy at two-body cluster level 




+ 



+ [ / d^r(7/;^(r)+i(;^(r))(l-e-^(^-^))(3(o-i-r) (0-2 • r) - o^i • 0-2) 



and show that 




Fourier transform of 512 (f) 



The term in © can be written in the form: 



(X yi2(0)-yi2(ki-kfc) , 



(D.l) 



where q = — kj and: 
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D. Tensor operator in momentum space 



The quantity in the first line vanishes, since the round brackets contain a function of 
|r| only. The one surviving term is the exchange contribution, which is proportional 
to the FT of the tensor operator. The integration is performed fixing the z-axis and 
denoting x — cos 6 = z - x' = cos 6' = z • q and q - f — cos 7 = y with q= |q|. 
We expand the exponential in Legendre polynomials according to 



21+1 /•! 
2 7-1 



(D.2) 



and use the standard decomposition of a scalar product v • r in spherical harmonics, 
as well as the addition theorem for Legendre polynomials 

1 

(Ti ' r 



U=-l 



47r 



p=-e 

where cr~^ denotes the spin rising and the lowering operator: 



1 




a- = V2 





1 



It follows that 

Vi2= j drr^W5{r)^ge{r,q) j dClPe{y)) 



1 



@ 



X 



Using the above property of Pe{y), the first term can be rewritten 

9e{r, q) ^ 



® 



= J drrV(r)Ef^(4-)^E(-)^^^'-r-2-''' = 
(I) 



^ dnYi^{x,<l>)Yi,,ix,<f>)Yi;ix,<l>) 
and, substituting 

Yi;{x,<f>) = {-rYi.p{x,cl,), 
we finally obtain the expression 

"9(2/ + 1)1^/2 / 1 1 / 

p 



(I) = E(-r 

p 

= E 



Aw 







1 1 I 

u u' —i 



9 



47r(2Z + 1) 



1/2 



(1100|/0) {llup'\lp) , 



7^ for I = 0, 2 

where the amplitude {Iil2mim2\lm) is the usual Clebsch-Gordan coefficient describing 
the addition of £i,mi and £2, "^2 in the £, m channel. Note that the first coefficient 
(1100|£0) is different from zero only for ^ = 0,2. 
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^ = term 

For i = the contribution denoted (a) turns out to be 



= drr^W3{r)go{r,q){47rf/^ 



5^00 (^S 00 



47r(cri-o-2) drr^w^{r)go{r, q) 



(Ti - (72 



1 

\/47r 



As a consequence, — © = in V12 and 

£ = 2 term 

In this case 

47r\3/2 



V^i2= /drr2(...)52(r,Q)(y) 3 ^(-r+''V^2"'^' (1100|20)( 1 1 i^i.' Kp) l2p(a;', <^'), 

uy'p 



or 



In order to simphfy the above result, we start from the definition 



and use 



{n00\m){lliyiy'\£p)Y2p 



In order to ehminate the i = contribution we add 



CTi • (72, 



to obtain the final result 



3(cri • q)(cr2 • q) - cri • cr2 



/6(47r) 



(_)-+-Vr%-^Vi^^^|2^)y2p, 
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D. Tensor operator in momentum space 



that has been employed in the calculation of the energy from 

- 5^12(4)^ J drr'^g2{r,q)(w^{r) +WQ{r)^ > 
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